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The short term variabihty of the Galactic cosmic ray flux (CRF) reaching Earth has been previ- 
ously associated with variations in the global low altitude cloud cover. This CRF variability arises 
from changes in the solar wind strength. However, cosmic ray variability also arises intrinsically 
from variable activity of and motion through the Milky Way. Thus, if indeed the CRF climate 
connection is real, the increased CRF witnessed while crossing the spiral arms could be responsible 
for a larger global cloud cover and a reduced temperature, thereby facilitating the occurrences of 
ice ages. This picture has been recently shown to be supported by various data In particular, 
the variable CRF recorded in Iron meteorites appears to vary synchronously with the appearance 
ice ages. 

Here we expand upon the original treatment with a more thorough analysis and more supporting 
evidence. In particular, we discuss the cosmic ray diffusion model which considers the motion of the 
Galactic spiral arms. We also elaborate on the structure and dynamics of the Milky Way's spiral 
arms. In particular, we bring forth new argumentation using HI observations which imply that 
the galactic spiral arm pattern speed appears to be that which fits the glaciation period and the 
cosmic-ray flux record extracted from Iron meteorites. In addition, we show that apparent peaks 
in the star formation rate history, as deduced by several authors, coincides with particularly icy 
epochs, while the long period of 1 to 2 Gyr before present, during which no glaciations are known 
to have occurred, coincides with a significant paucity in the past star formation rate. 
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I. INTRODUCTION 

It has long been known that solar variabihty is affecting 
climate on Earth. The first indication for a solar-climate 
connection can be attributed to William Herschel [Q, 
who found that the price of grain in England inversely 
correlated with the sunspot number. He later suggested 
that it was due to changes in the solar irradiance p5| . 
The irradiance variability is probably not large enough 
to explain the climatic variability observed by Herschel, 
nevertheless, synchronous temperature and solar varia- 
tions do exist. For example, typical surface temperatures 
during northern summers were found to differ by 0.5°K 
to 1.5°K between solar minima and solar maxima p5| . 

Over the past century. Earth has experienced a grad- 
ual, though non-monotonic warming. It is generally be- 
lieved to be a result of a greenhouse effect by anthro- 
pogenic fossil fuel emissions. However, a much better fit 
is obtained if part of the warming is attributed to a pro- 
cess, or processes, correlated with the solar activity, thus 
explaining for example, the non-monotonic global tem- 
perature change 1^, Moreover, the part of the 
climatic variability which is synchronized to the solar ac- 
tivity is larger than could be expected from just the 0.1% 
typical change in the solar irradiance H, ^ . Namely, the 
variability in the thermal flux itself appears to be insuf- 



* Electronic address: 



ihaviv@phys.huji.ac.il 



ficient to explain, for example, the global temperature 
variations observed. 

If one goes further back in time, then climatic variabil- 
ity on the time scale centuries, is too correlated with solar 
activity. Cold episodes in Europe such as the Maunder, 
Sporer and Wolf Minima clearly correlate with peaks in 
the flux, while warm episodes, such as the "medieval 
warm" period during which Vikings ventured across the 
Atlantic, correlate with minima in the ^''C flux (e.g., 
p7t). This flux itself is anti-correlated with the solar 
activity through the solar wind which more effectively 
reduces the Galactic cosmic ray flux that reaches Earth 
(and produces ^"^C) while the sun is more active. On a 
somewhat longer time scale, it was even found that cli- 
matic changes in the Yucatan correlate with the solar 
activity (and possibly with the demise of the Maya 
civilization). While on even longer time scales, it was 
shown that the monsoonal rainfall in Oman has an im- 
pressive correlation with the solar activity, as portrayed 
by the ^'^C production history [fz^ . 

Two possible path ways through which the solar activ- 
ity could be amplifled and affect the climate were sug- 
gested. First, solar variations in UV (and beyond) are 
non-thermal in origin and have a much larger relative 
variability than that of the total energy output. Thus, 
any effect in the atmosphere which is sensitive to those 
wavelengths, will be sensitive to the solar activity. UV 
is absorbed at the top part of the atmosphere (at typical 
altitudes of 50 km), and is therefore responsible for the 
temperature inversion in the Stratosphere. Any change in 
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the UV heating could have effects that propagate down- 
ward. In fact, there is evidence that it can be affecting 
global circulations and therefore also climate at lower al- 
titudes. For example, it could be affecting the latitudinal 
extent of the Hadley circulation . 

Another suggested proxy for a solar-climate connec- 
tion, is through the solar wind modulation of the galac- 
tic cosmic ray flux, as first suggested by Ney j?^. Ney 
pointed out that cosmic-rays (CRs) are the primary 
source for ionization in the Troposphere, which in return 
could be affecting the climate. 

First evidence in support was introduced by Tinslcy & 
Deen in the form of a correlation between Forbush 
events and a reduction in the Vorticity Area Index in 
winter months. Forbush events are marked with a sud- 
den reduction in the CRF and a gradual increase over a 
typically 10 day period. Similarly, Pudovkin & Verete- 
nenko reported a cloud cover decrease (in latitudes 
of 60N-64N, where it was measured) synchronized with 
the Forbush decreases. Later, an effect of the Forbush 
decreases on rainfall has also been claimed - an av- 
erage 30% drop in rainfall in the initial day of a Forbush 
event (statistically significant to 3a) was observed in 47 
Forbush events recorded during 36 years in 50 meteoro- 
logical stations in Brazil. While in Antarctica, Egorova 
et al. found that on the first day after a Forbush 
event, the temperature in Vostok station dramatically 
increased by an average of 10°K, but there was no mea- 
surable signal in sync with solar proton events. On the 
longer time scale of the 11-year solar cycle, an impressive 
correlation was found between the CRF reaching Earth 
and the average global low altitude cloud cover ^ . 

Although the above results are empirical in nature, 
there are several reasons to believe why the cosmic-ray 
route could indeed be responsible for a connection be- 
tween the solar variability and cloud cover. First, CRs 
are modulated by the solar activity. On average, the he- 
liosphere filters out 90% of the Galactic CRs (e.g., [[zsl ). 
At solar maximum, this efficiency increases as the solar 
wind is stronger. Since it takes time for the structure of 
the heliosphere to propagate outward to the heliopause 
at 50 to 100 AU and for the CRs to diffuse inward, the 
CR signal reaching Earth lags behind all the different in- 
dices that describe the solar activity (e.g., the sunspot 
number or the 10.7 cm microwave flux which is known 
to correlate with the EUV flux). The cloud cover signal 
is found to lag as well behind the solar activity and it 
nicely follows the lagging CRF. Second, both a more de- 
tailed analysis |^ and an independent study |7^ show 
that the correlation is only with the low altitude cloud 
cover (LACC). Among the different possible causes which 
can mediate between the solar variability and climate on 
Earth, only Galactic CRs can affect directly the lower 
parts of the atmosphere. It is the Troposphere where the 
high energy CRs and their showers are stopped, and are 
responsible for the ionization. EUV variability will affect 
(and ionize) the atmosphere at higher altitudes (> 100 
km). As mentioned, thermal heating by Ozone absorp- 



tion could possibly affect also low altitudes, however, the 
CRF-cloud cover connection is seen only in low altitude 
clouds. Solar CRs (which are less energetic than Galactic 
CRs) are not only stopped at similarly high altitudes, the 
terrestrial magnetic field also funnels them towards the 
poles. On the other hand, the LACC-CRF correlation is 
seen globally. 

Last, Tinsley & Deen who first found a correla- 
tion between Forbush events and a reduction in the Vor- 
ticity Area Index in winter months, showed that these 
events correlate significantly better with the cosmic ray 
fiux than with the UV variations (which generally start 
a week before the Forbush events). They also suggested 
that the UV cannot be responsible for this Tropospheric 
phenomenon since the time scale for the Stratosphere to 
affect the Troposphere is longer than the Forbush-VAI 
correlation time scale. 

Although the process of how CRs could affect the cli- 
mate is not yet fully understood, it is very likely that the 
net ionization of the lower atmosphere (which is known 
to be governed by the CRF, e.g., [Q) plays a major 
role, as the ionization of the aerosols could be required 
for the condensation of cloud droplets ||l8|, An 
experiment is currently being planned to study the pos- 
sible cosmic-ray flux - cloud-cover connection. It could 
shed more light and perhaps solidify this connection . 
Moreover, som e physical understanding appears to be 
emerging [115|. Interestingly, the latter work may ex- 
plain why the apparent effect is primarily on the lower 
troposphere and why the global warming of the past cen- 
tury has been more pronounced at the surface than at 
higher altitudes. 

Hence, the evidence shows it to be reasonable that so- 
lar activity modulates the cosmic ray flux and that this 
can subsequently affect the global cloud cover and with 
it the climate. Assuming this connection to be true, we 
should expect climatic effects also from intrinsic varia- 
tions in the CRF reaching the solar system. 

With the possible exception of extremely high energies, 
CRs are believed to originate from supernova (SN) rem- 
nants (e.g., Q, 1^). This is also supported with direct 
observational evidence ||2^. Furthermore, since the pre- 
dominant types of supernovae in spiral galaxies like our 
own, are those which originate from the death of massive 
stars (namely, SNe of types other than la), they should 
predominantly reside in spiral arms, where most massive 
stars are born and shortly thereafter die |2^. In fact, 
high contrasts in the non-thermal radio emission are ob- 
served between the spiral arms and the disks of external 
spiral galaxies. Assuming equipartition between the CR 
energy density and the magnetic field, a CR energy den- 
sity contrast can be inferred. In some cases, a lower limit 
of 5 can be placed for this ratio . 

Thus, when the sun passes through the Galactic spiral 
arms, an increased CRF is expected. If the CRF-LACC 
connection is real, this will increase the average LACC 
and reduce the average global temperature. The lower 
temperatures will then manifest themselves as episodes 
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during which ice ages can occur. Moreover, if the Milky 
Way as a whole is more active in forming stars, more 
massive stars will die and produce CRs. We show in this 
work that both these effects appear to be supported by 
various data. 

Shaviv p7| studied this conjecture and found evidence 
which supports it, thereby strengthening the possibility 
of a CRF-climate connection. In this work, we elaborate 
the original treatment by performing a more thorough 
analysis. We significantly extend the discussion on the 
dynamics of the spiral pattern of the Milky Way as it is 
important for determining the reoccurrences of ice age 
epochs, and introduce a new argument that helps deter- 
mine the pattern speed. We also introduce more evi- 
dence in the form of an apparent correlation between the 
recorded Milky Way activity (as described by the star for- 
mation rate, SFR) and the occurrence of ice age epochs. 
In particular, it is shown that the lack of glaciation ac- 
tivity on Earth between 1 and 2 Gyr BP (before present) 
appears to correlate with a dip in the Star formation rate 
in the same period (at least, as obtained by several but 
not all authors!). On the more speculative side, since the 
SFR activity may correlate with the activity in the LMC 
and with its estimated passages through perigalacticon, 
ice-ages could be attributed, to some extent, to fly-by's 
of the LMC. 

It is also interesting to note that other mechanisms 
have been previously proposed to link the Galactic en- 
vironment with climate variability on Earth. The first 
such mechanism was proposed long ago by Hoyle & Lyt- 
tleton who argued that an encounter of the Solar 
System with an interstellar cloud might trigger an ice 
age epoch by increasing the solar luminosity, which pro- 
duces an over compensating increase in cloudiness. The 
increased luminosity is a result of the accretion energy 
released. However, it is currently believed that radiation 
driving has a positive feedback (e.g., Q), not a strong 
negative one. Namely, an increase in the solar luminosity 
will result with an increase of the temperature, not a de- 
crease. Nonetheless, encounters with interstellar clouds 
could still have a temperature reducing effect if sufficient 
quantities of dust grains are injected to the uppe r atmo- 
sphere to partially shield the solar radiation |114|. These 
events are more likely to occur during spiral arm crossing, 
since it is there where dense molecular clouds concen- 
trate. However, since they require high density clouds, 
it seems unlikely that they can explain several 10^ yr 
glaciation epochs each spiral crossing. 

A second mechanism has to do with the shrinking of 
the heliosphere. Begelman & Rees Q] have shown that 
while crossing moderately dense ISM clouds with den- 
sities of 10^ to 10^ cm~^, the bow shock of the helio- 
sphere will be pushed further in than 1 AU. As a conse- 
quence, the slowing down effect that the heliosphere has 
on Galactic cosmic rays, will cease to work and the flux of 
Galactic low energy CRs will be significantly increased. 
On the other hand, the charged particles comprising the 
solar wind will not reach Earth. Either way, the flux 



of low energy charged particles reaching Earth could be 
significantly altered. Although these particles are not 
known to have a climatic effect at the moment, such an 
effect cannot be ruled out. Unlike the previous mecha- 
nisms, if this route can work, it may require significantly 
less dense ISM clouds which are more frequent. However, 
it is still unclear whether a several 10^ yr long glaciation 
event can be obtained via this route. 

A third mechanism operating mainly during spiral arm 
crossing is the perturbation of the Oort cloud and in- 
jection of comets into the inner solar system. Napier & 
Clube 1?^ and Alvarez et al. [Q discussed the effects that 
grains injected into the atmosphere by cometary bom- 
bardment will have on the climate by blocking the solar 
radiation. Hoyle & Wickramasinghc |Q proposed that 
a cometary disintegration in the vicinity of Earth's orbit 
would similarly inject grains into the atmosphere. 

One should note that these mechanisms all predict ice- 
age epochs in synchronization with the spiral arm cross- 
ing. This is counter to the model described here in which 
a phase lag exists. 

There were also proposals that related the Galactic 
year (i.e., the revolution period around the galaxy) to 
climate on Earth (p^,|[ll ^, [pO[ and references therein). 
For example, Williams |113|| suggested that lAEs on 
Earth are periodic, and that this rough ~ 150 Myr pe- 
riod is half the Galactic year. Williams raised the pos- 
sibility that this Galactic-climate connection could arise 
if the disk is tidally warped (e.g., by the LMC), but did 
not mention a specific mechanism that can translate the 
warp into a climatic effect. On the other hand, Steiner 
& Grillmair |^ suggested that climatic variability may 
arise if the solar orbit around the galaxy is eccentric and 
if, for some unknown physical reason, the solar luminos- 
ity is sensitive to the galactic gravitational pull. 

Another interesting suggestion for an extraterrestrial 
trigger for the ice-age epochs, has been made by Dilke 
& Gough [^9); (see also ||l^), who showed that the so- 
lar core may be unstable to convective instability under 
the presence of chemical inhomogeneities induced by the 
nuclear burning. These authors have argued that both 
the time scales and luminosity variations involved could 
explain the occurrence of lAEs. 

In addition to the extraterrestrial factors, there are also 
terrestrial factors which are in fact most often claimed 
by the paleoclimatological community to affect climatic 
variability on geological time scales. These are the con- 
tinental geography, sea level, atmospheric composition, 
and volcanic, tectonic and even biological activity. It is 
likely that, at least to some extent, many of the afore- 
mentioned terrestrial and extraterrestrial factors affect 
the global climate. Therefore, one of the main questions 
still open in paleoclimatology is the relevant importance 
of each climatic factor. 

We begin by reviewing the observations and measure- 
ments. These include a summary of the glaciation epochs 
on Earth, the dynamics and star formation history of 
the Milky Way, and the CRF history as derived from 
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TABLE I: Acronyms 



Notation Definition 



BP 


Before Present 


CR 


Cosmic Ray 


cR 


Co- Rotation 


CRF 


Cosmic Ray Flux 


HI, HII 


Atomic, Ionized Hydrogen 


lAE 


Ice- Age Epoch 


LACC 


Low Altitude Cloud Cover 


LMC 


Large Magellanic Cloud 


MW 


Milky Way 


SFR 


Star Formation Rate 


SN 


Supernova 



Iron/Nickel meteorites. Some of these results are de- 
scribed here for the first time. We then proceed to de- 
scribe the model which relates the Galactic environment 
to climate on Earth though the variability in the CRF, 
assuming CRs do affect the climate, and follow with the 
predictions of the model. The model's backbone is the 
solution of the problem of CR diffusion while incorpo- 
rating that the CR sources reside primarily in the spiral 
arms, and adding the climatic effect that the CRs may 
have. Then, we continue with a comparison between the 
proposed theory and observations. We show that an ex- 
tensive set of tests employing currently available data 
points to the consistency of the theory. 



II. THE OBSERVATIONS AND 
MEASUREMENTS 

We begin by reviewing several seemingly unrelated 
topics: The evidence for climatic variability on Earth on 
a time scale of 10^ — 10^ yr, as portrayed by the occur- 
rence of ice ages, the dynamics of the Milky Way with its 
spiral structure in particular, as well as the data on CR 
exposure ages from Fe/Ni meteorites. Some of the obser- 
vational conclusions are essentially quoted "as is" while 
several results are obtained by analyzing previously pub- 
lished data. 



A. Earth's glaciation history 

During the course of Earth's history, the climate has 
been variable on all time scales ranging from years to 
eons. Since clear geological signatures are left from pe- 
riods when Earth was cold enough to have extensive 
glaciations, studying the occurrence of ice-ages is a good 
method to quantify long term climatic variability, though 
it is not the only way (for example, we could have studied 
the occurrence of "evaporates" left during warm periods). 
We therefore choose to look at the occurrence of ice-ages. 

Before we continue, we should point out that ice-ages 



on Earth appear on two time scales. Over 0(10^) yr, 
ice ages come and go. However, epochs during which ice 
ages can appear or epochs during which ice ages can be 
altogether absent, exist on time scales of 0(10'') yr. In 
the rest of our discussion, the term ice-age epoch (lAE) 
will correspond to these long epochs. Today, we are in 
the midst of a long lAE, though specifically in a mid- 
glaciation period between 0(10'') yr long ice-ages. 

Extensive summaries describing the lAEs experienced 
by Earth are found in Crowell jl^ and in Frakes et al. 
P0| . These mostly rely on geological evidence of ice ages 
for the occurrence of glaciations, but not only. The na- 
ture of the glaciations in the Phanerozoic (0 — 545 Myr 
BP) are to a large extent well understood. Partially it 
is because more recent data is more readily available 
and partially because dating layers with foss ils is eas- 
ier. Moreover, the analysis of Veizer et al. [103| who 
measured the tropical sea surface temperatures over the 
Phanerozoic serves as an independent analysis from those 
studying the occurrence of glaciations. As can be seen 
in figure [l^ and table ^ the different analyses are quite 
consistent with each other. 

The Neo-Proterozoic (1000 - 550 Myr BP) was prob- 
ably an intrinsically cooler period in Earth's history, 
and glaciation was more abundant than in the Phanero- 
zoic. However, two epochs stand out as particularly more 
glaciated. These are ||l^ the Marinoan and Varangian 
Glaciations (545 - 585, 590 - 640 Myr BP) and the Stur- 
tian Glaciations (700 - 750 Myr BP). In the former, the 
extent of glaciations was particularly impressive, with 
evidence of low latitude sea level glaciations, which trig- 
gered ideas such as the 'Snowball Earth' |Q. A third, 
earlier epoch around ca. 900 Myrs BP is still very ques- 
tionable, with some less firm indications pointing to it 
940 Myr according to Wilhams |113t , and ca. 900 
Myr according to others, [Q, [0). To be conservative, 
we will not take this epoch in our analysis (though it 
does nicely correlate with a spiral arm crossing if it ex- 
isted). Before 1000 Myr BP, there are no indications for 
any glaciations, expect for periods around 2.2 - 2.4 Gyr 
BP and 2.9-3.0 Gyr BP. The lack of glaciations could be 
attributed to a changed solar orbit within the Galaxy. 
However, since the probability for the solar system to 
abruptly change its Galactic orbit is very small, this 
change which occurred at 1 Gyr BP, is more likely to 
be attributed to intrinsic variations in the climate — for 
example, due to a slow reduction in greenhouse gases, or 
to variations in the MW's average SN rate. 

The paleoclimatological data of Crowell Frakes 
et al. and Veizer et al. |103| is summarized in table 
H, together with our adopted age for the mid point of 
the ice-age epochs and its error. Panel C,D and E in 
fig. ^ depict a graphical summary of the appearance of 
glaciations in the past 1 Gyr. 

As a big word of caution, one should note that the 
glaciation data does not come without its caveats. For 
example, unlike Frakes et al. |3^, Crowell believes that 
the data is insufficient to claim periodicity in the oc- 
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TABLE II: Ice Age Epochs from geological records (in Myr 
BP) 

Midpoint of lAEs^ Adopted Age Spiral Arm Total 



Crowell Frakes Veizer 


and Error^ 


Error"^ 


Error 


< 15 


^ 28 


- 30 


20± 10 


11 


15 


~ 155 


~ 144 ' 


- 180 


160 ± 10 


11 


15 


~ 325 


~ 293 ' 


- 310 


310 ± 20 


12 


24 


~ 450 


~ 440 ' 


- 450 


446 ± 10 


13 


17 


~ 595 


~ 588 




592 ± 15 


14 


21 


~ 730 


~ 765 




747 ± 20 


15 


25 


~ 900 


~ 940 




920 ± 20 


16 


26 



^ Crowell |l^] is the mid point of the epochs with glacial 
activity. Frakes et al. is the midpoint of the "ice house" 
periods, while Veizer et al. [ lOJ ] is the time with the coldest 
tropical sea temperatures, for which evidence exists in the 
past 550 Myr. 

^ The adopted age is the aver age of Crowell jl^, Frakes 
et al. js^ and Veizer et al. [ lOJ ] , except for the present lAE. 
Since it is ongoing, its midpoint is likely to be more recent. 
The adopted error is a rough estimate. It considers that the 
error in recent lAEs is smaller and that short lAEs are much 
easier to pin point in time. 
^ The "Spiral Arm Error" is the error arising from the 
epicyclic motion of the solar system. That is, it arises from 
the non circular motion that it can have around the galaxy. 
It effectively introduces a 'jitter' in the predicted location of 
the spiral arms (see appendix B). 



currence of lAEs. See § V A and fig. |1C 
summary of the caveats. 



for a detailed 



B. Spiral Structure and Dynamics of the Milky 
Way 

The exact pattern speed of the spiral arms of Milky 
Way and in fact the spiral structure itself, is still consid- 
ered an open question. This is primarily because of our 
internal vantage point inside the Milky Way. Since these 
will soon be required, we review the current status and 
analyze the data available. 



1. Milky Way Spiral Structure 

A review of the different measurements for the spiral 
structure, and in p articular the number of arms is given 
by Vallee ^ |lOg and by Elmegreen H. Vallee con- 
cluded that 4 arms are more favorable than two. Elm- 



green concluded that 4 arms appear to govern the outer 
part of the Milky Way, while the inner part is much more 
complicated. The problem in the determination of the 
actual spiral structure is that distances to objects are ei- 
ther not known accurately enough or they do not trace 
the spiral arms unambiguously, if particular objects are 



used (for example, HII regions, OB stars or Cepheids). 
If a smooth component is analyzed instead (such as the 
distribution of molecular gas) then the distance, which is 
inferred from velocity measurements and the Milky Way 
rotation curve, is a multi-valued function of the gas veloc- 
ity within the solar circle. It is then hard to disentangle 
the spiral structure from the observed I — v (longitude ve- 
locity) maps. The main exception to the above is the HI 
(or similar) measurements of gas outside the solar circle. 
Since HI traces the spiral arms nicely, and since outside 
the solar circle no velocity-distance ambiguities exist, the 
spiral structure can be "read off" the I — v maps straight 
forwardly. The result is a clear 4-arm spiral structure 
pr] , |T^. We therefore assume that at the solar galac- 
tocentric radius and beyond, the spiral structure of the 
Milky Way is that of 4 arms. This does not imply that 
further inside the Galaxy the same 4-arm structure ex- 
ists. In fact, we shall show that there is a good reason 
for the two structures to be different (which could explain 
why until now the picture was confusing). 



2. The Spiral Arm Pattern Speed - Previous Results 

Even less agreed upon and more confusing are the re- 
sults for the pattern speed f2p of the spiral arms in the 
Milky Way[119|. For the sake of completeness, we first re- 
view previous determinations of Qp. We will afterwards 
continue with a new analysis which was previously ap- 
plied to other galaxies but never to our own. 

A survey of the literature reveals that quite a few dif- 
ferent analyses were performed to measure the Galactic 
spiral arm pattern speed. Some methods are local in 
the sense that they look at local age gradients of young 
objects, such as OB stars or open clusters. These are 
presumably methods which rely on the least number of 
(Galactic) assumptions. For example, they should detect 
the correct pattern speed irrespective of whether the spi- 
ral arm is a density wave or just a star formation shock 
wave (without a density wave associated with it), or irre- 
spective of whether the MW has 2 or 4 spiral arms. Un- 
fortunately, these methods tend to be inaccurate because 
of local "dispersions" and inaccuracies in age determina- 
tions. 

A second type of methods looks at the birth place of 
objects not as young as before. These include, for exam- 
ple, the birth place of open clusters with a typical age 



of a few 10 yr. This can in principle help place more 
accurate constraints on the pattern speed. However, un- 
like the previous methods, it requires a model for the 
spiral arms including their number, their amplitude and 
pattern speed, all of which should be fitted for. In re- 
ality, one often assumes both that the number of arms 
and their amplitude are given within the context of the 
density wave theory. Then, several different guesses for 
are guessed and the best fit is chosen. 
A third type of methods relies on fitting the observed 
velocities of stars to a spiral density wave, and in partic- 
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ular the non-circular residue obtained after their circular 
component is removed. The advantage of this type of a 
measurement is that it does not rely on age determina- 
tions at all, since it relies on the "instantaneous" config- 
uration. However, the residual kinematics are sensitive 
to the rotation curve chosen as well as to the spiral wave 
parameters. 

A fourth type of methods relies upon the identification 
of resonance features expected to arise from the spiral 
density wave theory. For example, Gordon [ p5[ identifies 
the observed discontinuity in CO emission at about 4 kpc 
with the inner Lindblad resonance. 

A summary of the vario us d eterminations of flp in the 
literature is found in table III . The main result apparent 
from the table is that most of the values obtained for 
0„ cluster within two ranges (Qp, — r2„ ~ 9 — 13.5 km 



sec~^ kpc~^and — fip ~ 2.5 — 5 km sec~^ kpc~^), 
with a third range being either a tail for the second or a 
real "cluster" of results (for which ^Iq — f2p ~ (-1) — (- 
4) km sec^^ kpc^^). Interestingly, the division between 
the clusters is not a function of the method used. For 
example, Palous et al. have shown that two equally 
acceptable values are obtained from the same analysis. 
Clearly, the results in the literature are still not con- 
verged, but possible values and unaccepted ranges can 
be inferred. 



3. Pattern Speed from HI Observations 

As previously mentioned, relying on particular objects 
to identify the spiral structure has the disadvantage of 
distance inaccuracies and that these objects do not al- 
ways trace the spiral arms nicely enough. One the other 
hand, mapping of various gas components has the disad- 
vantage that the velocity-distance ambiguities can com- 
plicate the analysis significantly. Therefore, analyzing 
gas at Galactic radii larger than the sun {R > Rq) has a 
clear advantage, as it avoids the above ambiguities and 
uncertainties. 

Blitz et al. found that a four armed|120| spiral 
structure in HI extends all the way to i?Hi,out ~ 2Rq. 
(Specifically, they found i?Hi,out = 20 kpc when taking a 
rotation curve in which Rq = 10 kpc and Vq = 200 km/s. 
For more up to date rotation curves, the value of i?Hi.out 
is lower but still roughly twice our galactocentric radius) . 
Moreover, more upto date HI maps, as traced by CO re- 
veal the spiral arms outside the solar galactocentric ra- 
dius even more clearly ||l^, thus reinforcing the results 
of Blitz et al. ||l^. This observation on the external ra- 
dius of the galactic spiral arms can be proven useful to 
constrain the pattern speed, if this 4-arm structure is a 
spiral density wave. 

According to spiral density wave theory, 4-armed spiral 
density waves can only exist within the inner and outer 
4:1 Lindblad resonances (e.g., (l^). Otherwise, the waves 
become evanescent. The inner and outer 4:1 Lindblad 



resonances, i?/4:i and Rqaa, are defined through: 

ilQ[HM;l) — iZp — ilQ[U0i:l) H , 

(1) 

where il{R) and k{R) are respectively the rotational fre- 
quency and the epicyclic frequency at radius R. There- 
fore, the constraint i?Hi,out ^ Rqaa that the arms should 
terminate before or at the outer Lindblad resonance can 
be rewritten as 



,0 



/^(^HI,out) 



(2) 



To obtain the numerical value of the r.h.s., we need 
to know the rotation curve of the Milky Way. We use 
the summary given by Oiling & Merrifield | |75| , which in- 
cludes the range of currently acceptable rotation curves 
for which Rq ranges from 7.2 to 8.5 kpc. For each rota- 
tion curve, we recalculate the outer extent of the HI arms 
which the Blitz et al. ||ll| result corresponds to. We then 
calculate the location of the resonances and their con- 
straint on the pattern speed (eq. ||) . This is portrayed in 
figure |l|. The results are given in table IV. They imply 
that: 

VLp < 16.9±2.5 km sec^i kpc^i (3) 
^G-^p > 9.1 ± 2.4 km sec~i kpc"\ 

Again, it should be stressed that the results assume 
that the spiral arms are a density wave. However, if they 
are, then the limit is robust since 4 spiral arms are clearly 
observed to extend to about twice the solar galactocentric 
radius. Spiral density wave theory has thus far been the 
most successful theory to describe spiral features in ex- 
ternal galaxies therefore, it is only reasonable that it 
describes the spiral arms in our galaxy as well. Moreover, 
in alternative theories in which the spiral arms are for ex- 
ample shocks formed from stellar formation, one would 
not expect to see spiral arms beyond the stellar disk, 
which is "truncated" several kpc inwards from i?Hi out 

MM , 

This last point, combined with the fact that HI is seen 
beyond the outer extent of the arms (to R ~ 3Rq) im- 
plies that their outer limit is probably the actual 4:1 res- 
onance. In other words, the relations given by eq. |^ are 
probably not just limits but actual equalities. The rea- 
son is that there is otherwise no other physical reason to 
explain why the arms abruptly end where they are ob- 
served to do so. (They do not end because of lack of HI 
nor does the stellar population has anything to do with 
it). 

Further theoretical argumentation can strengthen the 
last point made. The co-rotation (cR) point, has often 
been linked to spiral arm brightness changes in external 
galaxies. The reason is that the spiral arm shocks are 
important at triggering star formation. However, near 
co-rotation, the shocks are very weak. This is often seen 
in external galaxies as an "edge" to the disk, external to 
which the surface brightness is much lower (e.g., [p5|). If 
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the 4-arm pattern has the aforementioned pattern speed, 
then the cR radius can be predicted. This radius (rela- 
tive to Rq) is given in the seventh column of table 
We find that Rcr — Rq ~ 6.5 ± 1.5 kpc. For comparison, 
observations of stellar distributions show that there is a 
sharp cutoff of stars at RcR — Rq — 5.5 to 6 kpc according 
to Robin et al. Ml or, Rcr — Rq = 6.5 ± 2 kpc according 
to Ruphy et al. [p4[ . These numbers are also in agreement 
with a sharp truncation of the CO mass surface density 
at Rcr — Rq = 5 ± 0.5 kpc Q. Namely, observations 
are consistent with the cR radius predicted using the cal- 
culated pattern speed, provided that the sharp cut-off is 
related to the CR radius. This is also in direct agreement 
of Ivanov |^ who finds a co-rotation radius at 11-14 kpc 
for Rq = 8 kpc. 

The last argument for why the constraint given by eq. ^ 
should be considered an equality and not a limit is the 
following. If flp is significantly lower than the limit given 
by eq. 0, then the inner extent of the 4-arms, as con- 
strained by the inner 4:1 Lindblad resonance, should be 
outside the solar circle. However, according to the Blitz 
et al. ||ll[ and Dame et al. ||l^ data, the 4-arms extend 
inward at least to the solar circle. On the other hand. 



table III shows that the 4 arms cannot extend much fur- 
ther in. Thus, not only should eq. |^ be considered a rough 
equality, the spiral structure inside the solar circle has to 
have different kinematics than the structure of the ex- 
ternal 4-arms. This could explain the large confusion in 
the spiral structure and pattern speeds. This could also 
explain why several authors find that the solar system 
is located near co-rotation. If we are near CR, then it 
would be the co-rotation radius of the inner spiral struc- 
ture. Nevertheless, this point is still far from having a 
satisfactory explanation. See §VA for more details. 



Considering now that the first range of results in table 
III appears to be consistent with the density wave theory 
and the observations of HI outside the solar circle, we 
average the results in this range to get a better estimate 
for the pattern speed. We find Qq — rip = 11.1±1 
km sec~^ kpc~^. This translates into a spiral crossing 
period of 134 ± 25 Myr on average (taking into account 
the results of appendix B). 



C. Star Formation History of the Milky Way 



in the Galaxy. In other words, the SFR on long time 
scales, as recorded in nearby stars, should record long 
term changes in the Milky Way SFR activity. This may 
arise for example, from a merger with a satellite or nearby 
passages of one. 

Scalo using the mass distribution of nearby stars, 
found SFR peaks at 0.3 Gyr and 2 Gyr before present 
(BP). Barry Q and a more elaborate and recent analysis 
by Rocha-Pinto et al. |^ (see also references therein), 
measured the SFR activity of the Milky Way using chro- 
mospheric ages of late type dwarfs. They found a dip 
between 1 and 2 Gyrs and a maximum at 2-2.5 Gyrs BP. 
As a word of caution, there are a few authors who find 
a SFR which in contradiction to the above. More detail 



can be found in the caveats section SVA 



These SFR peaks, if real, should also manifest them- 
selves in peaks in the cluster formation rate. To check 
this hypothesis, the validity of which could strengthen 
the idea that the SFR was not constant, we plot a his- 
togram of the ages of nearby open clusters. The data 
used is the catalog of Loktin et al. |6^. From the his- 
togram apparent in figure ^ two peaks are evident. One, 
which is statistically significant, coincides with the 0.3 
Gyr SFR event. The second peak at 0.6 Gyr, could be 
there, but it is not statistically significant. Thus, we can 
confirm the 300 Myrs event. Note that this cluster his- 
togram is not corrected for many systematic errors, such 
as the finite life time of the clusters or finite volume ef- 
fects. As a result, the secular trend in it is more likely to 
be purely artificial. The same cannot be said about the 
non-monotonic behavior of the peaks. 

One source for a variable SFR in both the LMC and 
MW could be the gravitational tides exerted during LMC 
perigalactica. According to the calculations of the peri- 
galacticon passages, these should have occurred within 
the intervals: 0.2-0.5 Gyr BP, 1.6-2.6 Gyr BP and 3.4- 
5.3 Gyr BP (with shorter passage intervals obtained by 
Gardiner et al. |3^, and the longer ones by Lin et al. [ pO[ ; 
the further back, the larger the discrepancy). Interest- 
ingly, the 0.3 Myr and 2 Gyr BP SFR events are clearly 
located in the middle of the possible LMC perigalacticon. 



In general, the intrinsic flux of cosmic rays reaching 
the outskirts of the solar system (and which we will soon 
require) is proportional to the star formation rate (SFR) 
in the solar system's vicinity. Although there is a lag of 
several million years between the birth and death of the 
massive stars which are ultimately responsible for cosmic 
ray acceleration, this lag is small when compared with the 
relevant time scales at question. In the "short term", i.e., 
on time scales of 10® yrs or less, this "Lagrangian" SFR 
should record passages in the Galactic spiral arms. On 
longer time scales, of order 10^ yrs or longer, mixing is 
efficient enough to homogenize the azimuthal distribution 



Apparently, the MW activity is also correlated with 
SFR activity in the LMC. At 2 Gyr BP, it appears that 
there was a significant increase in the SFR in the LMC. 
Photometric studies of the HR diagram by Gallagher 
et al. j32| show a prominent increase in the SFR some- 
what more than 2 Gyrs before present. An increase in 
the SFR 2-3 Gyrs before present was also found by Val- 
lenari et al. |101|. Dopita et al. analyzed planetary 
nebulae and found that the LMC metallicity incre ased 
by a factor of 2 about 2 Gyrs BP, and Westerlund ]10g| ] 
has shown that between 0.7 and 2 Gyrs BP, the LMC 
had a below average SFR. 
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TABLE III: The observational determinations of the Milky Way spiral arm pattern speed. 



First SIq Q,p Q,q - Qp Method / Notes 

Author (km/s)/kpc (km/s)/kpc (km/s)/kpc 



Yuan (^iyoyaj 


zo 


~ io.o 


lie; 


Arm dynamics fit to Lin & Shu 


Yuan (iybybj 


O K 
ZO 


~ 13.0 


lie; 
^ li.O 


Migration of young stars 


LxOraon i^iy/oj 


ZO 


11 c: i_ 1 r 
11.0 It 1.0 


lo.O =t 1.0 


CO discontinuity at 4 kpc is ILR. 


r cLiu Lib I ± y ( ( 1 


ZO 


lO.O 


^ 1 _L . O 


v..yiLioLei uii Lii place 


Ljrivnev (^lyooj 


ZO 


1 9 1 
IZ-ID 


Q 1 Q 


Cepheid birth place 


ivanov (^lyooj 


97 K 


1 9n 

ID-ZU 


/ .0-11.0 


Cluster age gradients 


Comeron (1991) 


25.9 


16 ± 5 


10 ± 5 


Kinematics of young stars 


Creze (1973) 


~ 24.5 


~ 22 


2.5 ± 1.5 


Kinematics of young stars 


Palous (1977) 


25 


~ 20 


~ 5 


Cluster birth place 


Nelson (1977) 


25 


~ 20 


- 5 


Spiral Shocks profile & 21cm line 


Mishurov (1979) 


25 


23.6 ±3.6 


1.4 ±3.6 


Kinematics of Giants & Cepheids 


Grivnev (1981) 


25 


21 - 23 


2-4 


Kinematics of HII regions 


Efremov (1983) 


25 


18-25 


0-7 


Cepheid age gradients 


Amaral (1997) 


23.3 


21 ± 1 


2.3 ± 1 


Cluster birth place 


Avedisova (1989) 


25.9 


26.8 ± 2 


-1 ± 2 


Age gradient in Sag-Car 


Mishurov (1999) 


26 ±2 


28.1 ± 2 


-2 ±3 


Cepheids kinematics 


Fernandez (2001) 


25.9 


30 ±7 


-4 ±7 


OB Kinematics 


Here (HI) 26.5 ± 1.2 


16.9 ± 2.5 


9.1 ± 2.4 


rout ~ ro4:i of 4 HI arms^ 


Here (CRF var.) 




-10.5 ± 0.8 


Stat ± l'5gyg3 


CRF variability in Fe meteorites 


Here (ice-ages) 




40.4 ±0.3^ 


*stat ± l-5s^s3 


Fit to ice-age occurrence* 



^ Some results have no quoted error, but they should typically be ±1 km sec~^ kpc~^. For example, Palous et al (1977) 
checked specific pattern speeds: 11, 13.5, 15, 17.5, 20, 21.5 and found that only 13.5 and 20 agree with cluster birth places. 
^ In principle, the 4 HI arms can terminate at rni.out < ?'04:i (in which case flp can be smaller and Qq — Qp larger than the 
quoted numbers), however, as explained in the text, there is evidence pointing to rui.out actually being roi-.i- 
The origin of the systematic error is from possible diffusion of the solar system (both radially and along its azimuthal 
trajectory), relative to an unperturbed orbit. This is explained in appendix B. The values include the expected correction 

(0.54 ± 1.5 km sec^^ kpc^^) due to the solar metallicity anomaly. 
* The agreement between the bottom three results form the basis for the spiral arms - ice age epochs connection. 



D. Cosmic Ray Flux history from Iron Meteorites 

When meteorites break off from their parent bodies, 
their newly formed surfaces are suddenly exposed to cos- 
mic rays, which interact with the meteorites through 
spallation. The spallation products can be stable nu- 
cleotides, which accumulate over time, or they can be 
unstable. In the latter case, their number increases but 
eventually reaches saturation on a time scale similar to 
their half life. The ratio therefore between the stable 
and unstable nucleotides can be used to calculate the in- 
tegrated CRF that the meteorite was exposed to from 
the time of break up to its burning in the atmosphere. 
Generally, it is assumed that the CRF is constant, in 
which case the integrated flux correspond to a given age 
through a linear relation | |5^ , . 

A twist on the above, however, takes place in measure- 
ments employing the ^^K/^^K ratio. Since the unstable 
isotope in the pair has a half life slightly longer than 1 
Gyr, it does not reach saturation. By comparing the age 
of meteorites using this method to methods which em- 
ploy unstable isotopes with a short half life (decay time 



< a few Myr), it was found that the CRF in the past 
several Myr has been higher by about 30% than its av- 
erage over 150 to 700 Myr BP In principle, if the 
measurements were accurate enough, the slight inconsis- 
tencies between the two types of methods, as a function 
of time, could have been translated into a CRF history. 
However, a simulation shows that except for the flux vari- 
ation over the past several Myrs, this method becomes 
unfeasible. This can be seen in figure ^j. 

To extract the CRF, another method should be used. 
If we look at figure ^, we see that points on the graph 
which are separated by equal real-time intervals, tend to 
bunch near the minima of the CRF signal, when plotted 
as a function of potassium exposure age. This statistical 
clustering effect can be used to extract the CRF, if a large 
sample of dated meteorites is used. 

As previously mentioned, CR exposure dating assumes 
that the fiux had been constant in history. However, if it 
is not, the assumed CR exposure time will not progress 
linearly with the real time. During epochs in which the 
CRF is low, large intervals of time will elapse with only a 
small increase in apparent CR exposure age. As a result, 
all the meteorites that broke off their parent body during 
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FIG. 1: The dynamical frequencies of the Milky Way relevant to its spiral structure, plotted as a function of galactocentric 
radius. Using a typical rotation curve (taken from [[7a]), one can calculate O and the epicyclic frequency k. Spiral density wave 
theory states that a four arm spiral cannot exist forDp > n + k/4. Since a four arm spiral is observed in HI to extend out to 
'"HI, out ~ 27?0, an upper limit on Q.p can be placed. This Q.p implies that the inner edge of the arm structure should be located 
roughly at our galactocentric radius with only two arms possibly penetrating further in (e.g. [^), if an m = 2 component exists 
with the same fip. 



TABLE IV: Implications of the 4-arm HI^ data using spiral density wave theory and a possible range for the Milky Way rotation 
curve. 



Rq Vq 


Qq 




Rin 


Rin - Rq 


^cR — Rq 


^p,max 


— ^p,max 


b,min 


kpc ^ 


km 

s kpc 


kpc 


kpc 


kpc 


kpc 


km 
s kpc 


km 
s kpc 


km 
s kpc 


7.1 184 


25.92 


16.2 


8.3 


1.2 


5.1 


14.39 


11.52 


30.1 


7.1 200 


28.17 


15.7 


7.1 


0.0 


5.0 


17.64 


10.59 


29.5 


8.5 220 


25.88 


18.7 


6.4 


-2.1 


7.7 


18.21 


7.7 


27.4 


8.5 240 


28.24 


18.3 


7.4 


-1.1 


7.9 


19.35 


6.7 


27.2 


7.8 212 


26.5 


17.2 


7.3 


-0.5 


6.5 


16.9 


9.1 


28.6 


±0.7 ±30 


±1.2 


±1.5 


±1 


±1.7 


±1.5 


±2.5 


±2.4 


±1.5 



^ HI data is from Blitz et al. ||ll|. Rotation curves are from Oiling & Merrifield [ |75| |. 
^ The MW appears to have a bar which ends at 3 to 4.5 kpc. Since bars typically end between the 14:1 resonance and cR, a 
lower limit on the bar's pattern speed Qb, can be placed (assuming the bar ends at 4.5 kpc and coinciding with 14:1 
resonance). Since fip,max < fib.min, at least two different pattern speeds exist in the MW. 
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Logio t [yr] 



FIG. 2: The history of the Star Formation Rate (SFR). The squares with error bars are the SFR calculated using chromospheric 
ages of nearby stars [^. These data are corrected for different selection biases and are binned into 0.4 Gyr bins. The line and 
hatched region describe a 1-2-1 average of the histogram of the ages of nearby open clusters using the Loktin et al. j6^] catalog, 
and the expected 1-a error bars. These data are not corrected for selection effects (namely, the upward trend with time is a 
selection effect, favorably selecting younger clusters). Since the clusters in the catalog are spread to cover two nearby spiral 
arms, the signal arising from the passage of spiral arms is smeared, such that the graph depicts a more global SFR activity 
(i.e., in our Galactic 'quadrant'). On longer time scales (1.5 Gyrs and more), the Galactic stirring is efficient enough for the 
data to reflect the SFR in the whole disk. The dashed histogram underneath is the same as the histogram above it, though 
only with clusters having a better age determination {w > 1.0, as defined in Loktin et al. jHsj). There is a clear minimum in 
the SFR between 1 and 2 Gyr BP, and there are two prominent peaks around 0.3 and 2.2 Gyr BP. Interestingly, the LMC 
perigalacticon should have occurred sometime between 0.2 and 0.5 Gyr BP in the last passage, and between 1.6 and 2.6 Gyr BP 
in the previous passage. This might explain the peaks in activity seen. This is corroborated with evidence of a very high SFR 
in the LMC about 2 Gyrs BP and a dip at 0.7 - 2 Gyr BP. Also depicted are the periods during which glaciations were seen on 
Earth: The late Archean (3 Gyr) and mid-Proterozoic (2.2-2.4 Gyr BP) which correlate with the previous LMC perigalacticon 
passage |^^, and the consequent SFR peak in the MW and LMC. The lack of glaciations in the interval 1-2 Gyr before 
present correlates with a clear minimum in activity in the MW (and LMC). Also, the particularly long Carboniferous-Permian 
glaciation, correlates with with the SFR peak at 300 Myr BP and the last LMC perigalacticon. The late Neo-Proterozoic ice 
ages correlate with a less clear SFR peak around 500-900 Myr BP. 



this interval, will cluster together. On the other hand, 
epochs with a higher CRF will have the opposite effect. 
The CR exposure clock will tick faster than the real clock 
such that meteorites that broke off during this period will 
have a wide range of exposure ages. Thus, the number 
of meteorites per "apparent" unit time will be lower in 
this case. 



We use the data of Voshage & Feldmann [105 and 
Voshage et al. [106]. Together, we have a sample of 80 
meteorites which were '*^K/'**'K dated. In principle, we 



could also use meteorites which where exposure dated us- 
ing unstable isotopes other than ^°K, however, this will 
introduce more complications. Because the other iso- 
topes used for exposure dating have a relatively short life 
time, the relation between the real age and the exposure 
age is sensitive to the ratio between the recent CRF (over 
a few Myr) and the average over the past 1 Gyr. This is 
not the case with ^"K. The second advantage of using the 
^"K data is that it encompasses more meteorites, which 
is important for the statistics. In the future, a more ex- 
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tensive and detailed analysis should clearly consider the 
additional data and the complications that it introduces. 

Since we are assuming that meteoritic surfaces are 
formed homogeneously in time, we should be careful not 
to be fooled from the effects that real clustering can have. 
To minimize this, we removed all meteorites which have 
the same classification and are separated by less than 100 
Myr. These are then replaced by their average. 

The result is a data set containing 50 meteorites. We 
then plot a histogram of the CR exposure age of the 
meteorites in figure |4[ Immediately apparent from the 
figure is a periodicity of 143 ± 10 Myr. Moreover, if we 
plot a scatter plot of the error as a function of age, we see 
a tendency to have a higher error in points falling in the 
gaps between the clusters. This can be expected if the 
clustering signal is real since points with a smaller error 
will more easily avoid the gaps, thereby generating a bias. 
This consistency check helps assure we are looking at a 
real signal. If the apparent periodic signal would have 
been a random fluctuation, there would have been no 
reason to have larger errors in the random gaps, except 
from yet another unrelated random fluctuation. 

Next, we fold that data over the apparent period to 
overcome the systematic error arising from the slowly 
changing selection effects. Once we fold the data over 
the period, we find an even clearer signal. We perform 
a Kolmogorov-Smirnov test on the folded data and find 
that the probability for a smooth distribution to gener- 
ate a signal as non uniform as the one obtained, is only 
1.2% in a random set of realizations. Therefore, it in 
unlikely that the meteoritic distribution was generated 
from a purely random process. Conversely, however, the 
K-S test can only say that the distribution is consistent 
with a periodic signal. 

Furthermore, we can also predict the d istribution us- 
ing the diffusion model described in § tIIB , In the graphs 
we take the nominal case of D = 10^ cm^/sec and 
Ih = 2 kpc (as will be shown in §111), and find that 
the distribution obtained fits the predicted one. In par- 
ticular, the phase where the obtained distribution peaks 
agrees with observations. This implies that the signifi- 
cance is of order 1:500 (instead of 1:80) to have a ran- 
dom signal accidentally produce the apparently periodic 
signal and in addition has the correct phase. 

The agreement between the amplitude predicted and 
observed is less significant. The reason is that by chang- 
ing the diffusion model parameters, the contrast can be 
changed (as is apparent in table 0). 

Another point to consider is the error in the CR ages. 
Although errors were quoted wi th t he exposure age data, 
their values are quite ad hoc [104|. By comparing the 
Potassium age to ages determined using other methods, 
it is evident that the quoted errors over estimate the ac- 
tual statistical error. The scatter in the difference be- 
tween the Potassium age and the age in other methods is 
typically 30 Myr. This implies that the actual statistical 
error in the Potassium age determinations is at most 30 
Myrs, though it should be smaller since some of the error 



should be attributable to the second age determination. 
Since the error will have the tendency to smear the dis- 
tribution, the obtained contrast between the minimum 
and maximum flux is only a lower limit. 

To summarize, the Iron/Nickel meteoritic exposure age 
distribution appears to have a periodicity of 143±10 Myr. 
The signal is not likely to arise from a random process, 
and it also has the correct phase to be explained by spiral 
arm crossings. The CRF contrast obtained (in the folded 
distribution) is (max(<I>)/ min(<I>)) > 3. In addition it has 
been previously obtained that the flux today is higher by 
about 30% than the average flux over the past billion 
years. Both these constraints should be satisfied by any 
model for the CR diffusion in the galaxy. 



III. THE MODEL 

We now proceed with a detailed description for how 
the dynamics of the Milky Way is expected to affect cli- 
mate on Earth. We begin with the large astronomical 
scale and work our way down the chain of physical links. 
Namely, we begin with a model for the CR behavior in the 
Galaxy and its relation to the spiral arms in particular. 
This will enable us to quantitatively predict the expected 
variability in the CRF. We then proceed to estimate the 
effects that this variability will have on cloud cover and 
the effects that the latter will have on the average global 
temperature. 



A. Dynamics of the Milky Way 

The key ingredient which will be shown to be the "driv- 
ing" force behind climatic variability on the long ~ 10^ yr 
time scale, is our motion around the Galaxy. Variability 
arises because the Milky Way, like other spiral galaxies 
does not have cylindrical symmetry — it is broken with 
the presence of spiral arms. 

The basic parameters that determine this variability 
are the solar system's distance from the center of the 
Galaxy Rq , its Galactic angular velocity Qq , the angular 
pattern speed of the spiral arms Qp, and their number m. 
Given these parameters, a spiral crossing is expected to 
occur on average every time interval of 



At 



2tt 



spiral 



m\Qr. 



(4) 



For typical values (as obtained in §||) , a spiral crossing 
is expected to occur every 0(10^ yi's). The next step, is 
to understand and predict the CRF variability that will 
arise from these spiral arm crossing events. 
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FIG. 3; Theoretical comparison between different exposure ages of Iron meteorites and their real age. Plotted as a function 
of the Potassium exposure age (cik) are the real age (Oreai, in solid line) and a non-Potassium exposure age (aBe, such as using 
i0gy2i]sg-g (^g^^jjjg^ with a dashed line), and ok (using a dotted line, with a unit slope). Also plotted are the predicted CRF 
relative to the present flux {$/^o) as a function of Orcai, and pt-the (unnormalized) expected number of Potassium exposure 
ages per unit time, as a function of ok. A histogram of ax should be proportional to pt. The horizontal and vertical dotted 
lines describe how pt is related to the relation between Orcai and ok — equally spaced intervals in real time are translated into 
variable intervals in uk, thereby forming clusters or gaps in ok. The graph of ase vs. demonstrates that comparing the 
different exposure ages is useful to extract recent flux changes (which determine the slope of the graph). On the other hand, 
the graph of pt demonstrates that a histogram of aK is useful to extract the cyclic variations in the CRF, but not for secular 
or recent ones. 



B. A Diffusion Model for the long term Cosmic 
Ray Variability 



Qualitative theoretical and observational arguments 
were used to show that the CR density should be con- 
centrated in the Galactic spiral arms. The next step is 
to construct an actual model with which we can quan- 
titatively estimate the variability expected as the solar 
system orbits the Galaxy. 

The simplest picture for describing the CR content of 
the galaxy is the "leaky box" , which assumes no spatial 
structure in the CR distribution "inside" the galaxy, and 
that the probability for a GR particle to remain inside the 
galaxy falls exponentially in time, with a time constant 
Te.ib- Namely, the general equation describing the density 



of a specie i is given by (e.g., |^): 

dTY; _ ^ _^ ^ ^ / 1 ^ 1 ^ 1 \ 

dt ' ^ ' \Te,lb Ts Tr J 

^Q.-^. (5) 

Ti,cff 

Tj^off is the effective time scale for decay for specie i which 
includes loss from the galaxy, spallation destruction and 
radioactive decay. Qi is the sum of all sources for specie 
i which include its actual formation Qi and the result of 
spallation of more massive species j (which has a spalla- 
tion time scale of Tj and a branching ratio Pij to form 
specie i). 

Eq. H can be generalized to the case when spatial and 
energy homogeneity are not assumed. One then obtains: 

^ = DV^N, + A [b^E)N,] - ^ + g.. (6) 
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FIG. 4: Histogram of the Iron meteorites' exposure ages. The lowest panel marks the ok ages on the x-axis and the quoted 
age error on the y-ajds. Even by eye, the ages appear to cluster periodically. The second panel is a 1:2:1 averaged histogram of 
meteorites with a quoted age determination error smaller than 100 Myr, showing more clearly the clustering peaks. Altogether, 
there are 6 peaks from 210 to 930 Myr BP. The period that best fits the data is 143 ± 10 Myr. The third panel is similar to 
the first one, with the exception that the data is folded over the periodicity found. It therefore emphasizes the periodicity. A 
Kolmogorov-Smirnov test shows that a homogeneous distribution could generate such a non-homogeneous distribution in only 
1.2% of a sample of random realization. Namely, the signal appears to be real. This is further supported with the behavior of 
the exposure age errors, which supply an additional consistency check. If the distribution is intrinsically inhomogeneous, the 
points that fill in the gaps should on average have a larger measurement error (as it is 'easier' for these points to wonder into 
those gaps accidentally, thus forming a bias). This effect is portrayed by the dotted line in the panel, which plots the average 
error as a function of phase — as expected, the points within the trough have a larger error on average. 



The b{E) term corresponds to the slowing down of the CR 
rays. For the energies at interest, we can safely neglect it 
as its time scale is very long. Here, is the effective time 
scale for spallation or radioactive decay. The leaky box 
Te^ib should not be included explicitly in r^, since proper 
losses from the galaxy are implicitly included through 
losses from the boundary conditions. 

The simplest of such models are ID diffusion models 
which assume a slab geometry for the galaxy. More com- 
plicated models take more careful consideration of the 
structure of the galaxy by including the radial behavior, 
making them 2D in nature. However, these models do 
not consider that the sources reside primarily in the spi- 
ral arms (cf. Q). We therefore construct a simple model 
in which the spirals are taken into account. 

A cartoon describing the geometry of the problem 
solved in given in figure ^ In particular, we assume the 
following assumptions: 

1. The Galaxy is a slab of width 21 h- Within it, a 
diffusion coefficient D exists for the cosmic rays at 
the relevant energy ~ 10 GeV. Out of this region, 
the diffusivity is much larger such that CRs can 
effectively escape in a negligible time. This can be 
described with a boundary condition of the form 
^CR = at z = ±Ih- 



2. The CR sources are located in cylinders with a 
Gaussian cross-section. This deserves some expla- 
nation. First, in the Taylor & Cordes ||9^ model for 
the free electron distribution (which we use), the 
best fit to the free electron density is obtained with 
Cx = hz = 300 pc where ax is the width of the hor- 
izontal Gaussian distribution oc exp(— x^/2(T^) and 
hz is the width in the function sech^(z//iz) which 
was used for the fit. Namely, we simplify the prob- 
lem by assuming that the vertical distribution is 
Gaussian as well, but with the same length scale as 
in Taylor & Cordes 0. Second, the assumption 
that the spiral arms are straight cylinders is permis- 
sible, since the spirals are tightly wound: i « 12° 
with i being their pitch angle). That is, on the 
typical distance separating the arms, the radius of 
curvature of the arms is much larger. The actual 
distance taken between them is: 

dij = A0iji?0 sini, (7) 

where A(j)ij is the angular separation between arms 
i and j according to the Taylor and Cordes model. 
The distribution of the free electrons and of the 
cosmic ray sources should be similar since the mas- 
sive stars that undergo supernovae and produce 
CRs are also the stars that ionize the interstellar 
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medium. More careful consideration later on will 
actually show that there is a small but noticeable 
temporal lag between the two. 

3. The arms are moving at a speed vj_ — (ilp — 
f20)i?osini in a direction perpendicular to their 
axes. This is permissible since under the above ge- 
ometry, there is symmetry along the arms' axes, 
implying that their motion along their axis has no 
effect on the CR distribution. 

We should consider that we currently reside in the 
"Orion Arm" , which is only a spur or armlet . Al- 
though it is not part of the global structure of our galaxy, 
we are required to take it into account in the calculation 
of the recent CRF variation. As a consequence, we should 
expect because of it to witness a cosmic ray flux that is 
higher than predicted in the minimal four arm model. 
Since the density of HIT regions in this spur is roughly 
half of the density in the real nearby arms we as- 
sume it to have half the typical CR sources as the main 
arms. 

Since the lifetime of these "spurs" is expected to be of 
order the spiral arm crossing period (e.g., |^), they are 
not expected to repeat themselves after a whole revolu- 
tion, nor can we predict other possible "transient" spiral 
armlets that we might have crossed in the past. 

The free parameters therefore left in the model are 
the diffusion coefRcient D, the halo half width 1^, and 
the pattern speed of the spiral arms flp. Typical val- 
ues obtained in diffusion models for the CRs yield D ~ 
10^8 cmVs and i?/, ~ 2 kpc |, |l[ [08|. 

We assume that the heliospheric structure does not 
vary over intervals of a few 100 Myr. Even if it does, 
since the cosmic ray energies in question are large (> 10 
GeV/nucleon), the heliosphere will not affect the CRF 
by more than ~ 10% [Q. Thus, the change in the flux 
reaching the solar system is also the change in the flux 
reaching the Earth. 

We first solve for the CR flux itself and later solve 
explicitly for the Be nuclei which are used to date the 
"age" of the CRs. For the CR flux, we can neglect the 
effect of spallation. We also neglect the ionization energy 
losses. We therefore solve a simplified version of eq. ^ 



'dt 



= DV^* + S. 



(8) 



S is the source 
drical arms: 



E 



Ai exp 



distribution given by the moving cylin- 

where Ai and cr^.z are the amplitude of arm i (with a 
global normalization that is yet to be determined) and 
the horizontal and vertical width of the arms, x^^o — 
(j)i^oRQ sinz is the current day "location" of these arms in 
the model geometry in which the spirals are straightened. 
(j>i^o is the galactocentric angle of the arms at the solar 



galactocentric radius. These numbers, together with A(j)i 
are taken from Taylor & Cordes [0 . 

To solve the diffusion problem, we take full advantage 
of the fact that eq. || is linear. The problem is time 
dependent. Hence, we begin by writing the solution in 
the form: 

$(x,<)=/" dt'ifii^- {x,,o + v±t)±,t;t'), (10) 

^ — oo 

where (p{x,t;t') is the solution (at x, t) for the diffusion 
for the cosmic rays that where emitted at time t' . This 
function can be written as: 

(p{x, t;t') = ^ (^i(x - (xifi + v±t)Sc,t; t',i) (11) 



namely, as a contribution from each separate arm. We 
can then further use the method of mirror images which 
ensures that on the boundaries at ±Ih the CR flux will 
vanish. We do so by writing: 



m— — oo,oo 



(12) 

where 4>ifi is the flux at x and t of the CRs emitted at t' 
by arm i assuming it is centered at x' = at t' , when no 
boundaries are present. This has the form: 



= A, 



A, 



oo exp 
dx'dz' 

-oo 



(x — x')^ + (z — z 



4DAt 



— ) exp 



2ai 



inDAt 

^■^P \ 2{2DAt+<7l) ~ 2(2DAt+ai) ) 

27r^2L>At + al^2DAt + al 



2'KaxOz 



(13) 



with At = t — t' . From equation |T^, we have that on the 
z — plane (i.e., for \z\ <^Ih ^ 0{1 kpc)) 



ipi{x,z = 0,t;t',i) 



(14) 



e'^P (- 2(2I5ST^) ^4 (0, exp [-2ll/{2DAt + a^)] ) 



211^20 At + al^2DAt + al 



exp 



1 



( 2{2DAt+al)) 

^2T:{2DAt + al) Ih 

2.4674 / , 
DAt 



X exp 



'■H \ 



2-K{2DAt + (j1) 

r2 \ 1 



where is the Elliptic Theta function of the fourth kind 
and the approximation is accurate to ±3%. We find that 
the effect of the boundary condition is to introduce an 
exponential decay. This shows that the diffusion model 
with boundaries is very similar to the leaky box model, 
with a decay time constant of Te = If^/ {2A7D). In fact, 
for t — t' + ij1/{2D) > re/3 the error between a simple 



15 



exponential fit and the Elliptic theta function is less than 
10%. 

If we go back steps |l^ and we finally obtain: 

^ J-oo 27r^2D{t - t') + al 
^ ^4 (0,exp [-2ll/{2D{t-t') + al)\) 
^2D{t - t') + al 

(15) 

The global normalization of the amplitudes Ai is deter- 
mined by calculating the time average {Fcr{x, y = 0, t)) 
which from record in Iron Meteorites should be 28% less 
than todays flux 

1. CR lag after spiral arm crossing 

An interesting point worth particular note is that the 
CR distribution is expected to lag behind the spiral arm 
passages. Two separate physical mechanisms are respon- 
sible for such a lag. 

CR distribution skewness: Because the CR distribu- 
tion is skewed towards later times, the CRF is higher at 
a given time after the spiral arm crossing than before 
the crossing. This is because the CRs diffuse in the in- 
terstellar medium, but the spiral arms which are their 
source are moving, thus leaving behind them a wake of 
slowly diffusing CRs. Before an arm reaches the region 
of a given star, the CR density is low since no CRs were 
recently injected into that region and the sole flux is of 
CRs that manage to diffuse to the region from large dis- 
tances. After the spiral arm crosses the region, the CR 
density is larger since locally there was a recent injection 
of new CRs which only slowly disperse. 

This lag is intrinsic to the the diffusion model and 
therefore need not be considered separately. It can be 
seen in fig. ^ in the skewed distributions seen around 
each spiral arm crossing. Although the peak flux is lagged 
by only a small amount due to this effect, the skewness 
implies that the mid-point of epochs that are defined by 
the CRF beinglarger than a threshold value, will be 
lagged. Table ^ shows that this lag can range from 6 
to 19 Myrs, depending on the model (while considering 
only those that fit the observed ^"Bc/Be and CRF varia- 
tions) and assuming a threshold that will soon be shown 
to correspond to ice-age epochs. 

SN-HII lag: A second delay between the CRF distri- 
bution and the midpoint of the spiral arm passage arises 
because the definition of the spiral arm used does not 
have to coincide with the actual CR source distribution, 
which are SN remnants. In our case, the spiral arms lo- 
cation is defined through the free electron distribution as 
fitted for by Taylor & Cordes using the pulsar disper- 
sion data. The free electron density is primarily affected 
by luminous OB stars that are blue enough to ionize the 
interstellar medium. Therefore, the midpoint of the free 



electron distribution will occur around half the life span 
of the stars which dominate the ionization, namely, a few 
Myr since it is dominated by the most massive O stars 
that form. SNe on the other hand, occur at the end of the 
life of stars which can have lower masses and therefore 
have longer life spans than the stars that dominate the 
ionization. Since the least massive stars that undergo SN 
are about 8Mq and have a life span of typically 35 Myr, 
the average SN occurs at typically half this life span (as 
a result of the contribution of more massive stars). 

To better quantify this lag, we use the program of Lei- 
therer et al. [£9| which calculates various properties of a 
starburst population. In particular, it can calculate the 
flux of ionizing photons and the SN rate, both of which 
we require. We use standard parameters for a starburst 
Namely, 10^ Mq are formed at t — with a Salpeter 
IMF (i.e., a slope of -2.35), with a lower cutoff of 1 Mq, 
an upper cutoff of 100 Mq and a minimum mass of 8 
Mq needed to trigger a SN explosion. It also calculates 
the wind loss using a theoretical estimate and using the 
Kurucz-Schmutz spectra. The results are depicted in fig- 
ure ^ It shows that the average ionizing radiation is 
emitted at 2.0 Myr after the starburst, while an average 
SN takes place at 17.4 Myr after the starburst, giving an 
average delay of 15.4 Myr. The single most important 
parameter to which this lag is sensitive to is the lower 
cutoff mass needed to form a SN. On the other extreme, 
the global normalization does not affect the lag at all, im- 
plying that the total mass of stars formed and the lower 
cutoff for star formation are both unimportant. The fig- 
ure also shows a more realistic distribution obtained not 
for an instantaneous starburst, but one that has a Gaus- 
sian distribution with at = 16.5 Myr. This gives a Gaus- 
sian distribution for the ionizing radiation expected for 
an arm of width aarm ~ 300 kpc (for the pattern speed 
that will later be shown to agree with the various data). 
This width is the value obtained by the Taylor & Cordes 
p7| model. The result for the SN rate in this case, is 
similar to a Gaussian distribution with a new variance of 
at ~ 20 Myr. 

If we wish to incorporate these two results to the CR 
diffusion, we need to correct the following: 

1. The spiral arm locations, as defined by the CR 
sources, are lagged by 15.4 Myr after the spiral 
arms as defined through HII. 

2. The width of the CR distribution is spatially (or 
temporally) larger by a factor of 19.5/16.4 — 1.20 
than the one obtained from the Taylor & Cordes 
[0 model. 



2. Effects of inter-arm SNe 

We have assumed thus far that SNe occur only within 
the spiral arms. There are however two SN sources out- 
side the spiral arms. These are from infrequent 'field' OB 
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FIG. 5: The components of the diffusion model constructed to estimate the Cosmic Ray flux variation. We assume for simplicity 
that the CR sources reside in gaussian cross-sectioned spiral arms and that these are cylinders to first approximation. This is 
permissible since the pitch angle i of the spirals is small. 



stars and from SNe of type la, of which the progenitors 
are not massive stars. 

"Field" SNe: Most star formation in spiral galaxies oc- 
curs inside the spiral arms, as a result of gas being shock 
excited by spiral arm passages. Nevertheless, some star 
formation occurs outside the arms. To estimate the frac- 
tion of the latter, we look at two types of observations 
- the distribution of giant HII regions and of O stars. 
HII regions are formed from the ionizing radiation of the 
massive stars that later explode as SNe. According to 
Evans j2^, about 15% of these regions reside in the intcr- 
spiral regions. This can give an upper estimate for the 
inter-arm SNe rate (of types others than la), since 15% 
is the number fraction of regions — it is not weighed by 
the actual number of OB stars in each complex (which 
is expected to be higher in the arms, where HII clouds 
are larger on average). A lower limit for the SNe rate in 
the inter-arm regions would be given by the fraction of 
O stars that reside in the inter-arm regions. O stars are 
more massive of the stars that undergo SNe explosions; 
they are therefore expected to be more concentrated in 
the spiral arms. According to the data in Lynds [ p6[ , 
about 20 of the 400 or so nearby O stars reside outside 
the nearby spiral arms, or 5%. A reasonable estimate 
for the field SNe rate would therefore be 10% of the spi- 
ral arm rate. It is interesting to note that the discrep- 
ancy between SNe rate estimates from historical records 
and from observations in external galaxies has been at- 
tributed to SN being distributed mainly in the spirals of 
the galaxy |2l| . 

SNe Type la: SNe Type la are believed not to origi- 
nate from the core collapse of a single massive star. As 
a result, they can be found everywhere in the disk, not 
limited to within the spiral arms. For the same reason. 



they are also found in Elliptical galaxies. Their fractional 
rate in Spiral galaxies, depends on the spiral type and to 
a larger extent, the analysis used to derive the numbers. 
Typical numbers obtained are 23% for Sab and Sb's, and 
9% for Sbc-Sd ||, or, 8%-23% for Sab and Sb's, and 
10%-19% for Sbc-Sd |lO|. i.e., A typical number ex- 
pected for the Milky way is 15%. 

Both the above inter-arm SN sources add up to typi- 
cally 25% of the rate in arms. (And it is unlikely to be 
more than 35% or less than 15% of the rate in the arms). 
Does this ratio correspond to an additional constant in 
the CR source that is 25% of the CRs generated in the 
arms? It would if CR acceleration efficiency was the same 
in spiral arms and outside of them. However, theoretical 
argumentation points to an efficiency proportional to the 
ambient density n Q . There are also consistent observa- 
tions to corroborate this. Radio emission from SNRs in 
M33, which is presumably from CRs synchrotron emis- 
sion, is seen only in SNRs located in apparently denser 
environments |22|] . This would imply that the inter-arm 
SNe will accelerate less efficiently the CRs, by the typical 
density ratio between arm and inter-arm regions. 

Thus, the average inter-arm SN to arm SN CR gener- 
ation ratio, is expected to be: 



S = 



(Winter — arm) ( Pint c r — ar ni ) (/la + /ll,lb/i 



field OBs 



(•fa 



/ll,lb(l — /field OBs) ' 

(16) 

where /la is the fraction of SN la's, /n.ib = 1 — /la is 
the fraction of SNe of other type, and /field obs is the 
fraction of massive stars that form outside of the Galac- 
tic spiral arms. According to Lo et al. |6^, the typical 
arm/inter-arm density contrast in H2, is at least 3 but 
can be as high as 15. While it is typically 6 in Ha [p7[. 
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FIG. 6: The ionizing radiation flux and SN rate as a function of time for a starburst of 10 at time t = 0, assuming 
standard parameters |59]. The heavy weight hues describe the total H ionizing radiation flux (dashed) and the SN rate (solid). 
Since SNe are dominated by less massive stars than those responsible for the ionization, and since SNe occur only at the end 
of the stellar life, the SNe distribution is lagging behind the ionization distribution by 15.4 Myr (the average ionizing radiation 
is emitted at 2.0 Myr after the starburst, while the average SN explodes 17.4 Myr after it). The light weight lines are the same 
as before only the instantaneous starburst is replaced with a Gaussian distribution with a = 16.5 Myr, which reproduces the 
HII arms seen by Taylor and Cordes. The lag remains the same, but the distribution for the SNR widens to 20 Myr. 



Thus, for the nominal values for the different fractions 
and a density contrast of 6, we obtain that S ^ 5%. For 
the extreme values for the different fractions and density 
contrast, it can range from 1% to 15%. This is the "back- 
ground" constant CRF that one has to add to the spiral 
arm source of CRs. That is, when calculating the flux, 
we have to add a constant to it: 

$total = $arm + S ($arm) • (17) 



3. Results for the CR diffusion 

After incorporating the effect of the inter-arm SNe and 
the lag generated by the SN delay, results for the CRF 
can be obtained. Given D, v (-^ flp) and Ih, we can 
calculate the variability in the CR flux. We can also 
calculate the "age" of the CRs as will be measured by the 
^"Be/^Be ratio. This is described in Appendix A. Sample 
variations of the CRF and the Be age as a function of time 



are given in figure 0. 

4- Constraints on the diffusion model 

The primary constraint used to place limits on CR dif- 
fusion models is the ^"Be/Be ratio observed. This ratio 
depends on the typical duration the CRs spend between 
their formation and their detection in the solar system, 
i.e., on the CR's typical age. This is because Be isotopes 
are basically formed as spallation products. ^Be is stable 
and therefore accumulates with time. On the other hand, 
^"Be has a decay time scale of several Myr. It therefore 
saturates. 

In appendix A, the survival fraction of ^'^Be/Be is cal- 
culated using the spiral arm diffusion model. This ratio 
can then be translated to an effective confinement age 
that a "leaky box" model is required to have for it to 
yield the same Be ratio. The results are summarized in 
table M. This results should be compared with the mea- 
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FIG. 7: The cosmic-ray flux variability and age as a function of time for D = lO'^* cm'^/s and Ih ~ 2 kpc. The solid line is the 
cosmic-ray flux, the dashed line is the age of the cosmic rays as measured using the Be isotope ratio. The shaded regions at the 
bottom depict the location, relative amplitude (i.e., it is not normalized) and width of the spiral arms as defined through the 
free electron density in the Taylor and Cordes (1993) model. The peaks in the flux are lagging behind the spiral arm crosses 
due to the SN-HII lag. Moreover, the flux distribution is skewed towards later times. 



TABLE V: Results for the spiral arm CR diffusion model^ 
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(")min 


(a)max 


AT 




kpc 


Myr 


*0 
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•S>0 
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[Myr] 


[Myr] 


[Myr] 


[°K] 


0.32 


0.50 


10.54 


0.05 


2 


33 


44.66 


0.65 


5.32 


6.54 


4.75 


16.89 


32.40 


0.10 


0.50 


33.74 


0.07 


1 


96 


27.92 


0.62 


12.06 


24.98 


13.59 


75.74 


26.95 


0.03 


0.50 


105.43 


0.13 


1 


68 


12.93 


0.66 


17.19 


40.39 


17.61 


114.69 


22.13 


0.01 


0.50 


337.38 


0.17 


1 


60 


9.18 


0.71 


18.17 


41.45 


14.73 


128.06 


20.29 


1.00 


1.00 


13.49 


0.22 


1 


76 


7.87 


0.73 


6.64 


9.14 


7.19 


15.02 


21.75 


0.32 


1.00 


42.17 


0.23 


1 


58 


6.97 


0.70 


8.62 


30.90 


22.23 


58.88 


19.25 


0.10 


1.00 


134.95 


0.22 


1 


43 


6.62 


0.69 


17.37 


49.52 


29.77 


101.00 


17.21 


0.03 


1.00 


421.73 


0.21 


1 


48 


6.99 


0.72 


18.09 


50.10 


22.91 


121.00 


17.98 


3.20 


2.00 


16.87 


0.56 


1 


27 


2.26 


0.82 


5.95 


12.08 


11.59 


15.20 


10.02 


1.00 


2.00 


53.98 


0.54 


1 


23 


2.28 


0.80 


9.20 


38.12 


33.80 


50.68 


9.74 


0.32 


2.00 


168.69 


0.43 


1 


22 


2.86 


0.76 


10.71 


56.75 


45.18 


83.92 


11.22 


0.10 


2.00 


539.81 


0.29 


1 


26 


4.36 


0.73 


18.80 


59.73 


37.74 


107.39 


13.79 


10.00 


4.00 


21.59 


0.90 


1 


16 


1.29 


0.97 


3.70 


15.76 


15.76 


17.95 


3.66 


3.20 


4.00 


67.48 


0.84 


1 


12 


1.32 


0.92 


5.54 


42.52 


42.45 


48.92 


3.85 


1.00 


4.00 


215.92 


0.69 


1 


08 


1.56 


0.84 


9.28 
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57.54 
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5.46 


0.32 


4.00 


674.76 


0.49 


1 


13 


2.30 


0.78 


11.79 


66.62 


54.16 


90.75 


9.03 



Note: ^ The results summarized in this table are a subset of the all the models calculated (and appearing in flg. 



sured ^°Be/Be ratio. The various observed values for 
this number once translated to an effective CR age in 
the leaky box model yield 24j^g^ Myr (see appendix A). 
Namely, values ranging from 18 to 36 Myr are reason- 
able, and should be recovered by any consistent diffusion 



model, including the present one. 

Unlike other diffusion models, the spiral arm diffusion 
model can have more constraints placed on it using the 
observed CRF variability. Since previous models are all 
in steady state, they cannot under any condition, satisfy 
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FIG. 8: Models which fit the observational constraints. The free parameters in the diffusion model are the Diffusion coefficient 
D and the half width of the Galactic halo (relevant for the diffusion). "O" denote models which fulfill the constraint on the 
ratio between the maximum flux and the minimum flux (once folded over the 143 Myr period). "V" denotes models which 
satisfy the constraint on the ratio between today's flux and the average over the past billion years. "A" denote models which 
satisfy the "age" constraint of the CRs, as measured by the '^"B/B ratio. 



these constraints. 

Constraints on the CRF variabiUty can be placed us- 
ing CR exposure ages of Iron meteorites. Comparison 
between ages measured using '*^K/^''K dating in which 
the unstable isotope ''^K has a very long half life (1.27 
Gyr) and other measurements employing a "short lived" 
isotope (decay time < few Myr), have been consistently 
inconsistent. The only viable explanation being that 
the CRF in the past several Myrs has been roughly 
30% higher than the average over the past billion years 
|l|, g |§. We take $o/ ($) = 0.7 ± 0.1. An even more 
interesting constraint can be placed using the actual CRF 
variation, as ex tract ed from the distribution of CR expo- 
sure ages (see §11 D). This gives max ($)/ min($) > 3 . 
For various reasons given in § II D , an upper limit cannot 
be placed. 

Table and fig. | summarize the results. Values con- 
sistent with the different constraints on the effective leaky 
box age, with the ratio between the flux today and the 
average flux, and with the minimum to maximum flux 
contrast are emphasized in the figure. Since models are 
required to satisfy all three observed constraints, we can 



place limits on the diffusion coefficient D and the halo 
half width In- These are lO^^ cmVs <D< 10^^ cmVs 
and 0.5 kpc ;$ ^ 2 kpc (with the values being corre- 
lated). 

If we return back to table 0, we see that the models 
satisfying the constraints have a maximum variation of 
the CRF of (5 = (max$ - min $)/$today ^ 1-0 to 1.9 
and a lag oq of 6 to 19 Myr. These numbers will be 
needed later on when estimating the global temperature 
variations and the lag between spiral crossings and the 
occurrence of lAEs. 



C. The Cosmic- Ray Flux Temperature relation 

Wc proceed now to calculate the last link in the model, 
which is the relation between the CRF variations and 
temperature changes on Earth. We first calculate the 
effect directly and assume it takes place through LACC 
variations, as was empirically found by Marsh & Svens- 
mark [6^ . We then continue with an independent calcu- 
lation of this relation using global warming. This indirect 
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method does not assume that the CRF-cUmate connec- 
tion is through LACC variations, and it therefore serves 
as a consistency check. 

1. AT/A{CRF) through cloud cover variations 

We first calculate AT/A{CRF) by breaking the re- 
lation into two parts, into A{LACC)/A{CRF) and 
AT/A{LACC). 

The A{LACC)/A{CRF) relation: Since the CRF- 
LACC effect is on low altitudes clouds (< 3.2 km, JSTf), 
the effect should arise from relatively high energy CRs 
(> 10 GeV/nucleon) which can reach equatorial lati- 
tudes, in agreement with observations showing a bet- 
ter CRF-LACC correlation near the equator Thus, 
when estimating the forcing that the CRs have on cloud 
cover, the relevant flux is that of CRs that can reach a 
low magnetic latitude observing station and that has a 
high energy cut-off. The flux measured in the Univer- 
sity of Chicago Neutron Monitor Stations in Haleakala, 
Hawaii and Huancayo, Peru is probably a fair measure- 
ment of the flux affecting the LACC. Both stations are 
at an altitude of about 3 km and relatively close to the 
magnetic equator (rigidity cutoff of 12.9 GeV). 

The relative change in the CRF for the period 1982- 
1987 at Haleakala and Huancayo is about 7% while 
the relative change in the LACC is about 6% [ |67| . 
Namely, to first approximation, there is roughly a lin- 
ear relation between the relevant CRF (i.e., the equato- 
rial or low altitude CRF) and the cloud cover. Based on 
these observations, we now assume for simplicity that a 
linear relation exists between the CRF and the absolute 
LACC (with a current global average of 28%) throughout 
a range of ~ - 60% in the LACC. 

This assumption of linearity is at this point unavoid- 
able since at present, we don't have any theory to explain 
the physical process involved nor do we have any lab mea- 
surements to quantify its nonlinear regime. Thus, when 
a temperature estimate is obtained, we should keep in 
mind that this assumption of linearity was assumed, and 
that in principle, we can expect any correction factor of 
order unity. Moreover, at this point, no one can promise 
us that the cosmic ray effect on climate does not saturate 
at small corrections, which would quench any interesting 
"galactic" effect. 

The AT/A{LACC) relation: A change of 1% in the 
global LACC corresponds to a net effective reduction of 

the solar radiative flux of AF© 0.6 W/m^ ||||. Thus, 

a change of ±28% in the LACC corresponds to roughly 
a =f17 W/m^ change in the radiation flux. Without any 
feedback (such as cloud formation), an effective globally 
averaged change of 1 W/m^ can be calculated to corre- 
spond to a global temperature change of AT = 0.30°K. 
This rate of change Aq = AT/ AFq assumes that an aver- 
age surface area receives 70% (because of the 30% albedo) 
of a quarter (because Earth is a sphere) of the solar con- 
stant of Fq = 1367 W/m2. 



More detailed models give rates A which range from 
A = 0.3 to 1.1°K/(W m-2) ||, with typical values be- 
ing more like 0.7 — 1 Namely, there are positive 
feedbacks that increase the zeroth order relation by a 
factor of A/Ao = 1 to 3.7. Therefore, an assumed max- 
imum change of AF = =f17 W/m^ will correspond to a 
change of AT = =f5°K for A/Aq = 1, or to AT = =f18°K 
for A/Aq ~ 3.7. We shall take a nominal value of 
A — 0.85°K/(W m~^). This value corresponds to a de- 
crease of 0.14°K for an increase of 1% in the CR flux. 

Note that the value of A is a function of the time scale 
at question. On shorter time scales, the large heat capac- 
ity of the oceans acts as a moderator and a lower value 
of A is expected. 



2. AT/A[CRF) measured using global warming 

A second method for estimating the climatic driving 
force of CRF variations, is to estimate directly the re- 
lation between CRF changes and measurements of the 
global temperature change. When we do so, we should 
bear in mind that greenhouse gases too contribute to 
the global temperature variation. To decouple the two, 
we look at a period in which the global temperature de- 
creased. Since such a decrease cannot be explained by 
greenhouse gases, we can be safe that we are not mea- 
suring their contribution. 

During the period from the early 50 's until the early 
70's, the solar activity (averaged over the solar cycle) 
declined from a maximum to a minimum. According 
to the CRF-cloud picture, this resulted with a weaken- 
ing of the solar wind which increased the Galactic CRs 
reaching Earth, increased the cloud cover and reduced 
the average global temperature. And indeed, during this 
period, the average land and marine temperature in the 
Northern hemisphere has dropped by about 0.15°K. A 
more detailed analyses ^, Q which decomposes temper- 
ature trends into solar effects, anthropogenic and residu- 
als shows that the component attributable to solar vari- 
ability is actually larger - a reduction of about 0.2°K, 
where an increase of about 0.05°K is a result of human 
activity. 

To relate to Cosmic Ray flux variations, we use the ion- 
chamber data and the neutron monitor data in Climax 
quoted in Svensmark |9^ . Together with the ratio of the 
fluxes in Climax and in Huancayo/Haleakala one can 
obtain that the 0.2°K drop correlated with an increase 
of about 1.5% in the CR flux at Huancayo/Haleakala 
that presumably is responsible for the cloud cover effect. 
Namely, a 1% increase in the CR flux is responsible for 
a 0.13°K drop in global temperature. This result is sur- 
prisingly close to the estimate (of -0.14°K per 1% change 
in the CRF) that we obtained using the nominal value of 
the cloud-climate forcing. This is probably a coincidence 
since realistically, we should expect an uncertainty of a 
factor of 2 or so. Moreover, the long term effect of CRs on 
the temperature is probably somewhat higher since the 
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data used includes marine air temperature. Since the 
oceans have a high heat capacity (and therefore moder- 
ate the temperature), the long term temperature effects 
of CRs is expected to be higher. That is to say, -0.13°K 
per 1% change in CRF is probably a lower estimate. 



IV. THEORETICAL PREDICTIONS 
VS. OBSERVATIONS 

We now proceed with the comparison between the 
model predictions as derived in §111, and the observa- 
tional results summarized (and derived) in 



spiral arm asymmetries 

To check the consistency, a second analysis was per- 
formed to find the probability that a random distri- 
bution of glaciation epochs could generate a result 
which is as low as previously obtained. To do so, lAEs 
where randomly chosen. To mimic the effect that nearby 
glaciations might appear as one epoch, we bunch to- 
gether glaciations that are separated by less than 60 Myrs 
(which is roughly the smallest separation between ob- 
served glaciations epochs). The fraction of random con- 
figurations that surpass the obtained for the best fit 
found before is of order 0.1%. (If glaciations are not 
bunched, the fraction is about 100 times smaller, while 
it is about 5 times larger if the criterion for bunching is 
a separation of 100 Myrs or less). 



A. Spiral Structure vs. Climate on Earth 



For the nominal values chosen in our diffusion model, 
the expected CRF varies from about 25% of the current 
day CRF to about 135%. This should correspond to an 
average absolute LACC of 7% to 38% respectively. For 
the nominal value of A/Aq = 2.8, it should correspond to 
a temperature change of +10°K to — 5°K, relative to to- 
day's temperature. This range is definitely enough to sig- 
nificantly help or hinder Earth from entering an ice-age. 
(The typical difference between today, an interglacial pe- 
riod and the recent ice ages is of order 5°K). The fact 
that we are currently in a glaciation epoch and that the 
observed CRF is higher than the long term average, over 
which Earth was mostly out of glaciation epochs, is con- 
sistent with this picture. 

Comparison between the cosmic ray flux and the oc- 
currence of glaciations shows a compelling correlation for 
the pattern speed chosen (fig. |l^). To quantify this cor- 
relation, we perform a analysis. For a given pattern 
speed of the arm, we predict the mid-point of the spiral 
arm crossings using the data of Taylor & Cordes 1 97 , as 
described by Leitch & Vasisht |Q . The mid-point of the 
glaciation epochs are then predicted to lag by 15.4 Myr 
(from the SN-HII lag) + 6-19 Myr (from CRF skewness) 
« 21 — 35 Myr after the mid-point of the spiral crossing. 
Fig. y describes the value of — f2p which minimizes 
the fit between the predicted and measured mid-point 
of the lAEs, and the minimum value of jv obtained 
(with = 6 is the number of degrees of freedom minus 
the number of free parameters). This is done as a func- 
tion of the assumed lag between the mid-points of the 
spiral crossing and the predicted lAE. We find that the 
minimum ^ jv of 0.5 is obtained for Tiag = 33 Myr. This 
is consistent with the predicted lag. On the other hand, 
the graph shows that the no lag case is rejected at the 2cr 
level (85%). This implies that lAEs appear to be lagged 
after spiral crossings as predicted. On the other hand, the 
graph shows that placing the arms symmetrically (with 
a 90° separation instead of the values obtained by Tay- 
lor & Cordes |9^) does not improve or degraded the fit. 
It only shifts the preferred lag to slightly larger values. 
Hence, this cannot be used to learn about the possible 



B. Comparison with CRF record in Fe/Ni 
meteorites 



A smoking gun, which would unequivocally demon- 
strate that the spiral-arm — > cosmic-ray climate con- 
nection is real, would be a "historic" record which would 
reveal that the CRF was indeed variable. Such an 
independent record could on one hand show that the 
CRF varied as expected from the spiral arm CR dif- 
fusion model, while on the other hand, that it varied 
synchronously with the appearance of ice-age epochs on 
Earth. In §IID, a new method was developed with which 



we could extract the actual CRF history. 

It was found that the CRF record shows a clear vari- 
ability signal with a periodicity of 143 ± 10 Myr. This 
variability is consistent both in periodicity and phase, 
with the predictions of the CR diffusion model which 
incorporates the spiral arms, of which the expected peri- 
odicity is 134 ± 25 Myr. This is clearly demonstrated in 

Since the periodicity in the CRF is more accurately 
known than the spiral arm pattern speed (because of both 
statistical and systematic errors), it is a much stronger 
tool to use when comparing the astronomical data with 
the geological one. Indeed, the periodicity of the lAEs as 
can be obtained from table H is found to be 145.5 ± 7 
Myr, which nicely agrees with the aforementioned CRF 
period of 143 ± 10 Myr. This agreement in phase can also 
be seen in fig. Q where the clusters of exposure ages tend 
to be anti-correlated with the occurrence of lAEs. 

The last result is of course trivial once we consider that 
the CRF as recorded in Fe/Ni meteorites was already 
shown to agree with the predicted CRF using the galactic 
model and that the latter was shown to agree both in 
phase and period with the occurrence of lAEs. 

The implications are not only statistical, but also qual- 
itative, since the above agreements provide a "missing 
link" which explicitly point to the cosmic rays as being 
the culprit in climatic variability. Namely, it looks as if 
we found the smoking gun we sought after. 
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arms and the middle of the glaciation epoch (solid line). The hatched region describes expected lag from theory (including 
both SNR-HII lag and CRF skewness). The short dashed line is — f2p (and its error, when in the shaded region) which best 
fits the glaciations as a function of the lag. The long dashed line is a repeat for x 1^ when assuming a symmetric arm location 



The past Star Formation rate history and 
ice-epochs 



A sum mary of the past star formation rate history was 
given in § p C| . It is summarized again in figure ^, together 
with the epochs during which Earth witnessed ice-ages. 

The main correlation apparent from the graph is the 
paucity in the MW's star formation rat e bet ween 1 and 
2 Gyr BP (with the cautionary note of §V A), which co- 
incides with a long interval on Earth during which there 
were no apparent glaciations whatsoever. On the other 
hand, during the last eon, and the one between 3 and 
2 Gyrs BP, Earth did have glaciations and the SFR was 
higher. Although it is speculative, this variability in SFR 
appears to correlate with perigalacticon passages of the 
LMC, and in fact, also with the SFR in the LMC. 

For the SFR on shorter time scales, we have to rely 
on the cluster data (see fig. ^), or the mass distribution 
of nearby stars fs^ , which have a 'finer' resolution than 
the SFR calculated using chromospheric ages of nearby 
stars. These show that at 0.3 Gyr BP, there was a higher 
than average SFR. A second peak appearing in the clus- 



ter data, could have occurred at 600 Myrs BP (though 
the data is not statistically significant). Interestingly 
enough, the lAE 300 Myrs and 600-700 Myr BP appear 
to be correlated with more extensive glaciations: The 
recent peak correlates with the Carbonaceous-Permian 
lAE, while the peak at 600-700 Myrs coincides with the 
late Neo-Proterozoic lAE, both of which were relatively 
extensive. 

Although it is hard to quantify the statistical signifi- 
cance of the SFR-climate correlation, the fact that the 
results are consistent with that expected from the t heor y 
is reassuring (though see the possible caveats in § V A ). 
Namely, we expect and indeed find the relation: Higher 
star formation rate — > more CR generation — s- colder cli- 
mate. 



V. DISCUSSION 

The basic evidence presented supports a picture in 
which climate on Earth is affected by our changing lo- 
cation in the Milky Way, by way of a variable cosmic-ray 
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FIG. 10: Earth's recent history. The top panel (A) describes past crossings of the Galactic spiral arms assuming a relative 
pattern speed of f2p — = —11.0 km sec"'^ kpc~'^(which best fits the lAEs). Note that the Norma arm's location is actually 
a logarithmic spiral extrapolation from its observations at somewhat smaller Galactic radii The second panel (B) 

describes the Galactic CRF reaching the solar system using the CR diffusion model, in units of the current day GRP. An 
important feature is that the flux distribution around each spiral arm is lagging behind spiral arm crossings. This can be 
seen with the hatched regions in the second panel, which qualitatively show when lAEs are predicted to occur if the critical 
CRF needed to trigger them is the average CRF. Two dashed lines mark the middle of the spiral crossing and to the expected 
mid-glaciation point. Panels (G) (D) and (E) describe the paleoclimatological record of the past Eon. The solid line in panel 
(G) d epic ts the tropical sea surface temperatures relative to today, as inferred from calcite and aragonite shells in the past 550 
Myr 1 103 1 . The flUed areas describe the paleolatitudinal distribution of ice rafted debris [103|. Panel (D) and (E) qualitatively 
describes the epochs during which Earth experienced ice-ages, the top part as described by Frakes et al. , while the bottom 
one by Growell ||l^. The Phanerozoic part of panel E is directly taken from Growell ||l^. The Proterozoic part is its extension 
using the data in Growell. It has an elevated base level (since glaciations were more common) and it describes the main epochs 
of extended glaciations. By flne-tuning the actual pattern speed of the arms (relative to our motion) to best fit the lAEs, 
a compelling correlation arises between the two. The correlation does not have to be absolute since additional factors may 
affect the climate (e.g., continental structure, atmospheric composition, etc.). Panel (F) is a 1-2-1 smoothed histogram of the 
exposure ages of Fe/Ni meteors. The meteor exposure ages are predicted to cluster around epochs with a lower GRF flux. 



flux. The various physical links and the observational evi- 
dence which support these relations are charted in fig. |l^. 
One should note the redundancy in the key results which 
support this conclusion. These are the temporal correla- 
tions between the expected spiral arm crossing, the vari- 
able CRF observed in Iron meteorites, and the appear- 
ance of ice-ages on Earth. 

For nominal values taken, we find an effect of -1-5 to 
-10°K relative to today, however, there are some large 
uncertainties. A factor of order unity arises from the un- 
certainty in the cloud - temperature relation. Another 
factor of order unity arises from the uncertainty in the 
relation between the Cosmic Ray flux itself and the cloud 
cover. The normalization of the latter can be achieved us- 
ing the solar modulation with the complication however 
that it modulates CR fluxes of different energies differ- 
ently. Nevertheless, it is clear that the low altitude nature 
of the CRF-climate effect, requires the relevant energies 
to be high (> 10 GeV). Thus, a normalization factor of 
2 higher or lower between the CRF and cloud cover con- 



nection is conceivable. A consistency check which verifies 
that the net CR LACC temperature relations that 
we obtained are reasonable, is using global warming in 
the past century to directly calculate the CR temper- 
ature relation. 

The last source for an inaccuracy in the CR — > tem- 
perature relation arises because we extrapolate a several 
percent effect to several tens of percent. Without know- 
ing the physical process behind the effects of the CRs on 
the atmosphere. We cannot know, for example, whether 
the effect does not saturate at small relative changes. 

Although the actual links presented here relate the 
CRF to the changing global temperature, there is good 
reason to believe that the link itself is through a vari- 
able cloud cover. First, such a link was already empir- 
ically observed in the form of a direct correlation be- 
tween the low altitude cloud cover and the CRF. More- 
over, the observed correlation yields a radiation driving 
that is more than sufficient to explain the occurrence of 
ice-ages, but also the changes in the past century which 
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are attributable to solar activity. In other words, other 
pathways relating the CRF with temperature are not re- 
quired at this point, but clearly, we cannot rule them 
out. 

If the link is not through the variable low altitude cloud 
cover (or if an additional link exists) , then lower energy 
CRs could be affecting the climate as well. In this case, 
the effects that a dense ISM cloud has on the heliosphere 
could clearly be affecting the climate as well since 
these events will drastically change the low energy charge 
particle flux reaching Earth. In such a case, the ice- 
age epochs could be dotted with relatively short events 
during which the climate changes dramatically for the 
typical durations it takes to cross such ISM clouds (either 
heating or cooling, depending on whether the reduced 
flux of solar CRs or the increased galactic CRs are more 
important, and on the actual climatic mechanism). 

For example, if relatively low molecular cloud densi- 
tites of 100 cm^ can significanty alter the heliospheric 
structure and with it the low energy cosmic ray flux 
reaching Earth, and if the latter do affect climate, then 
an effect could actually take place for 10'' years, which 
is the typical duration it would take the solar system to 
cross a 100 pc cloud. In any case. These events are more 
likely to be synchronized with the spiral arm crossing 
itself and not lagged behind it. 

One should also note that passing through young stel- 
lar nurseries could also locally increase the intrinsic CRF 
reaching the solar system. Through the low altitude 
cloud cover connection, this would cause dramatic short 
term reduction in the temperature. Namely, the long 
ice-age epochs could be dotted with short events of even 
colder climate. 

The CR diffusion model developed here, has shown 
that the spiral structure should be incorporated in CR 
diffusion models since qualitative and quantitative cor- 
rections are otherwise missed. (The model described here 
is nowhere as complex as some CR models, so neither it 
nor previous models should be the final word). Moreover, 
even the spiral structure itself deserves a more detailed 
description. For example, the spiral arms were assumed 
to be straight "rods". A much bigger effect could how- 
ever arise from the approximation that the diffusion co- 
efficient is the same everywhere. However it is expected 
that the diffusion coefficient will be lower in the denser 
spiral arms. This could have an effect on the entrap- 
ment of the CR in the spirals for a longer time, thereby 
increasing the CR density contrast. 

To the extent that spiral density wave theory applies to 
the Milky Way's spiral arms, the spiral structure of the 
outer parts of our galaxy (R > Rq) is apparently pinned 
down. Therefore, it should serve as an outer boundary 
when Milky Way models are constructed. We found that 
the structure between us and the center will however be 
different in nature. The reason is that the outer 4-arm 
spiral cannot extend much further inside the solar galac- 
tocentric orbit, though spiral arms are seen further inside 
the disk. This implies that models previously used to de- 



scribe the galaxy with only one pattern speed are not 
sufficient. It could explain why previous results for the 
number were in contradiction. 

Another interesting observation is that geology is now 
is a unique position to serve as a tool in the study of 
galactic dynamics. For example, it now appears that 
the life time of spiral arms is at least ~ 900 Myr or 1.5 
rotations. This is expected whether the spiral arms are 
a 'transient' phenomenon or a more permanent one (an 
issue which is still not settled) . Should however the same 
periodicity be found in the glaciations 2-3 Gyr BP, then 
a very interesting constraint could be placed on the spiral 
dynamics. 

Moreover, using temporal variability (cither in the 
CRF or the glaciations), we have managed to measure 
the spiral pattern speed more accurately than using pure 
'astronomical' methods. The accuracy with these tem- 
poral methods is in fact large enough that one needs to 
consider that the corrections from our diffusion in the 
galaxy results in an offsetted pattern speed, since this 
correction is of order the error bars themselves or even 
larger. This is the first time a temporal record (instead 
of a "snapshot" ) is used to study galactic dynamics. 

It is also worth noting that the effects of terrestrial 
magnetic reversals on climate are probably not dramatic 
as could be naively expected, even though the total CRF 
is significantly increased in such events. The main reason 
is that the LACC-CRF connection is observed in low al- 
titudes and low latitude clouds. This necessarily implies 
that high energy CRs are those mostly responsible for 
the effect. These CRs reach the equator irrespective of 
the strength of Earth's magnetic field. For example, if 
we look at the original results of Compton |l^, we find 
that the change in atmospheric ionization at sea level 
varies by only ~ 14% between the magnetic equator and 
poles. Therefore, even if the Earth's magnetic field were 
to completely vanish, our estimates for the temperature 
change would correspond to at most a 1 to 2°K reduction 
in temperature. This is corroborated by the Laschamp 
event (36-41.5 kyr BP) during which the terrestrial field 
was significantly reduced Jsst . During the event, an in- 
creased ^"Be and '^^Cl flux was recorded in Greenland ice 
cores, but no statisticall y sig nificant climatic effect was 
recorded in ^'^O or CH4 [107|. This is expected consider- 
ing that the temperature variability recorded in the ice 
cores, either before or after the event, is much larger than 
1 to 2°K. 



A. Caveats 

The picture presented here includes many elements 
which add up to what appears to be a coherent pic- 
ture. Unfortunately, most if not all observational data 
used here, whether astronomical, meteoritic, or paleocli- 
matological does not come without a grain of salt. The 
purpose of this section is to discuss and summarize the 
limitations of the data and possible caveats. Hopefully, 
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FIG. 11: The model for how the Galaxy affects the ice ages. Solid arrows reflect physical cause and effect. Bold lines are hard 
to dispute relations while thin lines are reasonable relations based on observational correlations. Arrows with question mark 
are unknown though possible relations. Dotted and dashed arrows reflect observational correlations which relate two (or more) 
physical ingredients. Dotted lines are previously found correlations while dashed lined are correlations obtained and described 
here. 



by illuminating where the problems may rest, future re- 
search will try to address these points with the goal of 
resolving them. 

Paleoclimatological Data: By far, the least conclu- 
sive data used is the past evidence for the occurrence of 
glacial epochs, fn particular, it is debated whether ice- 
age epochs in the past billion years have been peri odic or 



not. Some claim that a periodicity exists ||30| , 113|, while 
others claim that insufficient data exists to claim period- 
icity Besides the timing of the glaciations (which 
becomes poorer as we go back in time), there is also the 
problem that ice-age epoch severity clearly varies from 
one lAE to the next. For example, the mid-Mesozoic ice- 
age epoch (~ 150 — 170 Myr BP) was clearly very 'mild' 
compared to other lAEs. On the other hand, from 800 
Myr BP until well after 600 Myr BP, it seems that Earth 
continuously had some sort of glaciations present, with 
colder episodes around 600 and 750 Myrs BP. In fact, 
around 750 Myrs BP, there is evidence for glaciations at 
sea level and equatorial latitudes. That is, glaciations 
more extensive than anything seen by Earth at any pe- 
riod p8[ . In addition to the problem of a varying severity 
or extensiveness, there is also the question of length. For 
example, the Carboniferous-Permian glaciations were 3 



or 4 times the duration of the short Silurian-Devonian 
glaciations that preceded. These problems are also noted 
in figure |l^. Note however, that even though the above 
problems exist, the timing of the lAEs of Frakes and of 
Crowell are generally consistent with each other. 

Another point which one might note at first glance of 
fig. ^ is that although the general agreement is good, 
there are some inconsistencies between spiral crossings 
and actual occurrence of lAEs. As previously explained, 
this should not be a concern because of several reasons. 
First, the actual spiral crossing signal is expected to have 
a 'jitter' because of the epicyclic motion of the solar sys- 
tem. Second, since the Milky Way is not a 'grand design' 
spiral galaxy, the spiral arms are not the perfect geo- 
metrical objects they are assumed to be. In particular, 
'gaps' in which the spiral arms are weaker or arm 'spurs' 
which fill the inter-spiral-arm with young CR producing 
stars are not uncommon. The lifetime of the spurs, how- 
ever, is less than a galactic orbit. Therefore, they are 
not expected to repeat themselves after 4 spiral cross- 
ings. Third, geological factors are also expected to be 
important in climate change. For example, the changing 
continental structure might free the poles of land masses, 
and could even allow temperate ocean currents to reach 
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the poles. This can hinder Earth from triggering the for- 
mation of ice sheets. 

The above factors appear to be particularly important 
in the last lAE. Although the main spiral arm crossing 
took place about 50 Myr BP, the last lAE started only 
'recently', about 30 Myrs BP. The reason we are now 
in an lAE is easy to explain since we are in the Orion 
armlet. This leaves the question of why didn't the lAE 
begin as the solar system entered the Sagittarius-Carina 
arm? This could be a result of factors intrinsic to Earth. 
The main factors often cited as the reason that the latter 
part of the Cretaceous was one of the warmest periods in 
Earth's history (e.g., ||5^), are that the land mass did not 
occupy the poles and that the global sea level was highest 
than ever (thus moderating the temperatures). Neverthe- 
less, if careful attention is given to the temperature vari- 
ations over the Cretaceous and Cenozoic (e.g., |118 ) 
then it is apparent that the temperature variations were 
more complicated than a simple monotonic cooling. At 
90 Myr BP, the average global temperature peaked. But 
in the 10 Myr preceding the K/T transition (at 66 Myr 
BP), the temperature dropped by more than 5°K. The 
temperature rose again by several degrees, and peaked 
at around 52 Myr BP, after which it declined (though 
non-monotonically) until present, when extensive glacia- 
tions are common. It is tempting to attribute the first 
minimum to the crossing of the Sagittarius-Carina arm. 
Namely, Earth did experience a reduced global temper- 
ature however, because of intrinsic effects (and perhaps 
because the arm at the point of crossing was not promi- 
nent), the temperature did not fall enough to trigger 
glaciations. After the Sag. -Car. arm was crossed, the 
temperature rose again. However, the Orion arm and 
the changing global conditions caused a large cooling af- 
terwards. 

Although the above explanation is far for satisfactory, 
it shows that the model developed in the paper is incom- 
plete. In particular, it is probably the case that "stan- 
dard" paleoclimatic factors intrinsic to Earth are as im- 
portant at determining the global climate as the external 
effect of the cosmic rays. 

Astronomical Data: Wc have used in the analysis 
astronomical data describing the spiral structure, spiral 
pattern speed and star formation rate. Unfortunately, 
there is still no consensus on any one of them. 

Although the analysis presented here, and in particular 
the meteoritic exposure age data, seems to be consistent 
with the various spiral arm indicators located outside the 
solar galactocentric radius, if we look inwards, the struc- 
ture and dynamics appear significantly less clear. In par- 
ticular, there are convincing arguments which indicate 
that the solar galactic orbit is close to co-rotation (e.g., 
ifzol ). If true, it would imply a spiral arm pattern speed 
similar to Qq, instead of a much larger Hq — J7p, as found 
here. Obviously there is a contradiction that should be 
resolved. One possibility raised in §2, is that the Milky 
Way has more than one set of spiral arms with different 
pattern speeds. This is required if the outer parts have a 



low pattern speed, and it would explain the large confu- 
sion that exists in the field. Clearly, this issue is still not 
resolved and deserves more attention than it has received 
in recent years. 

Interestingly, having two sets of spiral arms also has 
the potential of explaining the 'staggered' severity and 
length seen in the ice-age epochs. Namely, it could be 
that we are also affected by the spiral structure internal 
to our galactocentric radius, and that its period rela- 
tive to us is about 300 Myr. Nevertheless, this is quite 
speculative, and it is not unlikely that this staggered be- 
havior arises from terrestrial factors. For example, at- 
mospheric CO2 levels around the Carboniferous-Permian 
period were several times smaller than the other periods 
in global history with the exception of recent geological 
times. 

Another caveat has to do the location of the Norma 
arm at our galactic radius. This location was obtained 
through a logarithmic extrapolation from where it is ob- 
served at smaller galactic radii (a la Leitch & Vasisht 
p8|, T aylor & Cordes |Q). However, as discussed in 
i pl B 3| , there are good reasons to believe that the struc- 
ture within the solar circle is different from the one at 
our radius and outwards. This would imply that the 
Norma arm seen at smaller radii is not necessarily part 
of the same structure that should be localized between 
the Scutum-Crux and the Perseus external arms, on the 
other side of the galaxy. If so, the preferred crossing of 
this arm would be later by ~ 40 Myr, as its preferred 
location would be symmetrically in between the adjacent 
arms. If we look at fig. |l^, we see the this will in fact 
increase the agreement between the predicted location 
of the spiral crossings around 300 and 900 Myr BP and 
both the time of occurrence of the relevant lAEs and the 
observed CRF variability. It would not explain though 
the extent the the Carboniferous-Permian lAE. 

A third astronomical caveat has to do with the star for- 
mation rate. The SFR with which the glacial activity on 
Earth was compared with was that of Rocha-Pinto et al. 
p3[ , which is also consistent with various other authors. 
However, like most results quoted here, there are a few 
contradictions published in the literature, which, if veri- 
fied will pose a problem for the model presented here. In 
particular, Hernandez et al. |Q have used the Hippar- 
COS astrometric data combined with detailed theoretical 
isochrones to obtain the SFR of nearby stars (which due 
to diffusion corresponds to the azimuthally averaged SFR 
on periods over a billion years). The gross features they 
obtained are inconsistent with other analyses (for exam- 
ple, with ^^). In particular, the SFR they obtained de- 
creased from 2 Gyr BP until today with only relatively 
small 'wiggles'. Moreover, they did not find a dip be- 
tween 1 and 2 Gyr BP. Thus, an independent study of 
the HiPPARCOS data should be pursued in order to re- 
solve this discrepancy. 

Meteoritic Data: 

Regarding the variable CRF from exposure age data, 
one should note the following. Just the data showing the 
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clustering of meteorites by itself can be interpreted in 
several ways. The clustering can either be a manifesta- 
tion of a variable CRF, or, it can be a genuine clustering 
of exposure ages. 

Since the stable isotopes used in the currently em- 
ployed methods for exposure dating Fe/Ni meteorites 
(e.g., ■^^Ar, ^^Ne, ''^K), are all formed primarily from 
Galactic CRs the observed data could not be a re- 
sult of a variable solar CRF. Thus, within the interpreta- 
tion of a variable CRF, the observed variability can arise 
from either an intrinsic variability in the Galactic CRF 
(as assumed in this work) or a variability induced by a 
variable heliosphere. In the second possibility, one can 
accredit the exposure data to a variable solar wind (with 
the observed periodicity of ^ 150 Myr) and attribute 
the observed correlated climatic variability to the solar 
variations. Although possible, this option which has far 
reaching implications, seems less reasonable because of 
the correlation of the CRF and climatic data with the 
observed galactic pattern speed with the same period. 
Namely, it would be strangely coincidental if the sun were 
to have a variability that is synchronized with the spiral 
arm crossings. 

The third possibility for interpreting the exposure age 
data is that the CRF was constant and that the clus- 
tering in the exposure ages is real. Namely, episodes in 
the past in which asteroids or comets were more likely to 
break apart existed periodically. This is not an unrea- 
sonable hypothesis. For example, one could claim that 
the Oort cloud is perturbed more often during spiral arm 
crossing, injecting more comets to the inner solar system, 
thus breaking apart more asteroids. Moreover the more 
numerous debris could reasonably have climatic effects 
as well (see §|). Within the context of this interpreta- 
tion, clusters of meteoritic exposure ages, as well as the 
ice age epochs, should coincide with the spiral arm cross- 
ing with no phase lag. The data however shows that the 
clusters are not in phase with either the spiral crossings 
or the lAEs, nor are the lAEs in phase with the spiral 
arm crossings. 

Thus, the most consistent picture arises if we interpret 
the exposure age clusters as a result of an intrinsically 
variable Galactic CRF. One way to prove this, is to look 
at age distributions of Iron meteorites which were CR 
exposure dated using a method in which the stable iso- 
tope that accumulates is sensitive primarily to the solar 
cosmic-ray flux. The reason is that if the clustering is real 
(and not due to a variable Galactic CRF) then the same 
clustering will appear in solar sensitive dating methods. 
This would also be useful to rule out the possibility that 
the apparent variations are not Galactic in origin, but 
are induced by solar variability instead. If this unlikely 
case is true (which would imply that lAEs arc linked 
with solar variability, since the CRF-climate correlation 
exists irrespectively), then the opposite clustering will be 
present in solar sensitive data. With the most reasonable 
interpretation, though, no clustering is predicted in the 
solar sensitive data, once it could be obtained. 



VI. SUMMARY AND CONCLUSIONS 

The thesis presented here relates the following top- 
ics: The Milky Way dynamics. Cosmic Ray diffusion in 
the Galaxy, the CR record in Iron meteorites, the effect 
of Cosmic Rays on climatology, and glaciology. This is 
achieved by considering the intimation that CRs can af- 
fect the global cloud cover and that the CR flux from the 
Galaxy should be variable. Before summarizing the main 
conclusions relating the above topics, "preliminary" con- 
clusions on each of the various topics can be drawn from 
the background analysis. 

Milky Way Dynamics: By studying the 4-arm struc- 
ture seen in HI, extending to roughly 2Rq , stringent con- 
straints can be placed on the spiral pattern and dynamics 
of the Milky Way: 

1. Since observations of HI outside the solar circle suf- 
fer no velocity- distance ambiguity, the observations 
of 4-spiral arms should be considered robust. // one 
further assumes that these arms are density waves 
(which is by far the most consistent explanation for 
spiral arms) then these arms cannot extend further 
out than the 4 to 1 Lindblad resonance. This im- 
plies rip < 16 ± 2.5 km sec~^ kpc~^. This number 
agrees with roughly half the determinations of Qp, 
which cluster around this value. 

2. Several additional considerations point to the above 
value being not an upper limit for Op, but the ac- 
tual value of it. 

3. The 4-arm spiral cannot extend much further in 
from the solar galactocentric radius. Several pos- 
sibilities were raised. Since the MW has a bar, at 
least two different pattern speeds exist in the MW. 

4. On a more speculative note, long term activity in 
the star formation rate appear to correlates with 
activity in the LMC and possibly with its peri- 
galacticon passages, suggesting that long term star 
formation in the Milky Way and LMC are related. 

CR Diffusion: 

5. Since CR sources are clearly more common to the 
Galactic spiral arms, models for CR propagation in 
the MW should take this fact into account. Oth- 
erwise, significant discrepancies can arise. In par- 
ticular, quantities such as the "age" of the cosmic 
rays reaching Earth can be distorted. 

6. For typical diffusion model parameters, the CRF is 
expected to vary by 0(1), which is consistent with 
radio observations of external spiral galaxies. 

7. The distribution of CRs is both lagging behind the 
spiral crossings (when defined for example by the 
HII regions) and skewed towards later times. The 
skewness arises from the asymmetry introduced by 
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the spiral arm motion. The lag in the peaks arises 
because the SN explosions are lagging the ionizing 
photons that produce the HII regions, and which 
were used here (and in ||9^) to define the location 
of the arms. 

CR record in Iron meteorites: 

8. The "standard" method for extracting historic vari- 
ability of the CRF, by comparing ^^K/'^'^K dating 
to dating with a short lived unstable nucleotide 
(such as ^°Be/^^Ne) are effective only at extract- 
ing "recent" changes (over several Myrs) or secular 
changes over longer durations, but they are not ef- 
fective at extracting the signal expected from the 
periodic variability of the spiral arm passages. 

9. An effective method for finding the periodic CRF 
variations is a statistical analysis of CR exposure 
ages. It assumes that the rate at which "new ex- 
posed surfaces" appear in Iron meteorites does not 
vary rapidly. 

10. This method shows that the exposure age data 
is consistent with it being periodic, with a pe- 
riod of 143 ± 10 Myr, and a CRF contrast of 
max(<I>)/ min(<I>) > 3. It uses 50 meteorites which 
were dated only with '^^K/^^K. The method reveals 
that statistically, it is unlikely that the meteor ex- 
posure ages were generated by a random process. 

The Milky Way Climate connection: By compar- 
ing the above findings to the appearance of glaciations 
on Earth, the following conclusions can be deduced: 

11. Periodic glaciation epochs on earth in the past 
1 Gyr can be consistently explained using the 
proposed scenario: Periodic passages through the 
Galactic spiral arms are responsible for an in- 
creased CRF, an increased LACC, a reduced global 
temperature and consequent ice-ages. 

12. Long term glaciation activity is apparently related 
to the global SFR activity in the Milky Way This 
may be related to LMC perigalacticon passages 
which would imply that the nearby extragalactic 
environment could be added to the factors affect- 
ing the global climate. 

The evidence upon which the above conclusions were 
based, is the following: 

I The period with which spiral arm passages have oc- 
curred using the HI data alone (and which should 
also correspond to the CRF period), is 163 ± 50 
Myr. By combining the result of the 8 total mea- 
surements that scatter around this value, the pre- 
dicted spiral arm crossing period is 134 ± 22 Myr. 
The period with which the CRF appears to vary us- 
ing Fe/Ni meteorites, is 143 ± 10 Myr. The period 



with which glaciations have been observed to reoc- 
cur on Earth is 145.5 ± 7 Myr on average. Clearly, 
all three signals are consistent with each other. 

II By comparing the actual prediction for the location 
of the spiral arms to the glaciations, a best fit of 
143 ± 5 Myr is obtained. The phase of all the three 
signals using this periodicity is then found to fit the 
predictions. In particular: 

• The average mid-point of the glaciations is 
lagging by 33 ± 20 Myr after the spiral arm 
crossing, as is predicted (31 ± 8 Myr, from a 
possible range of ±75 Myr). 

• The CR exposure ages of Iron meteorites clus- 
ter around troughs in the glaciations. 

III A random phenomenon for the appearance of 
glaciations is excluded with very high statistical 
significance. 

IV Long term variability in the appearance glaciations 
correlates with the observed SFR variability of the 
Milky Way. In particular, the lack of apparent 
glaciation between 1 and 2 Gyr BP, correlates with 
a particularly low SFR in the milky way (less than 
half of today), while a high SFR rate between 2 
and 3 Gyr BP (with a peak towards 2 Gyrs), cor- 
relate with the glaciations that Earth experienced 
between 2 and 3 Gyrs BP. 

Additional conclusions: More conclusions can be 
reached by considering the above results. They do not 
pertain to climatic variability, but instead, they can be 
drawn by assuming that long term climate variability is 
indeed a measure of CRF variability, as the above results 
seem to endorse. 

13. Once more accurate "Eulerian" measurements of 
the Galactic pattern speed will be available, inter- 
esting constraints could be placed on solar migra- 
tion in the galaxy. This is achieved by comparing 
the "Eulerian" measurements of the pattern speed 
to the "Lagrangian" pattern speed (as measured by 
the moving solar system), which is now known to 
an accuracy of ±3%. 

14. The variability observed in the CRF record in Iron 
meteorites (see point 10) can be used to place con- 
straints on models for the CR diffusion in the Milky 
Way. This is in addition to the currently used con- 
straint on the survival ratio of spallation products 
(primarily Be). Thus, from CRF variability, we find 
D < 2 X lO^s cmVsec and < 2 kpc. 

15. The 4-arm spiral pattern has been stable for at 
least a billion years. A better understanding of 
past Glaciations (e.g., 2-3 Gyrs ago) could place 
interesting constraints on the lifetime of the spiral 
arms which cannot be done by any other means. 
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As apparent from the above list, there are several inter- 
esting ramifications to the picture presented here. How- 
ever, one which bears particular interest on global warm- 
ing and which was not discussed here, is that it now ap- 
pears even more plausible that cosmic rays indeed affect 
the global climate, as suggested for example by Ney @, 
Tinsley & Deen 1 98 , and Svensmark & Friis-Christensen 
psf . This implies that it is now more plausible that so- 
lar variations are affecting climate through modulation 
of the cosmic ray flux. This would explain an impor- 
tant part of the global warming observed over the past 
century. It is because of this important aspect, in par- 
ticular, that we are obligated to iron out the still un- 
clear points raised in the caveats section, whether it be 
the paleoclimatological record, the possible mechanism 
by which cosmic rays affect climate, or our astronomi- 
cal understanding of the structure and dynamics of the 
Milky Way. 
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In the spiral arm diffusion model, there is no steady 
state, so the calculation of the ratio Nj/Ni is more com- 
plicated. 

To simplify the calculation, we utilize the previously 
obtained result that to a reasonable approximation, the 
effects of the boundaries is to have the CR distribution 
decay exponentially, not unlike the leaky box model. 
(Namely, we approximate the i?4-function as an expo- 
nent). This implies that the number density of the CRs 
that were emitted at a given time decays exponentially 
with a time constant of Tf. = 1\/{2A7D). Thus, the 
density of specie k which is the parent nucleotide of the 
spallation products (namely, C, N and O producing the 
Be) and which where generated at time t' will decay as 

exp(-(t-i')/Te). 

The two (united) equations that describe the produc- 
tion and decay of the spallation products formed by par- 
ent nucleotides that were generated at time i', are there- 
fore: 
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dt 
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(19) 



ni j{t, t') is the density of specie i (stable) or j (unstable) 
that were generated from nucleotides of type k aX t — t' . 

The solution to the two equations with the assumption 
that riijit = t') = (namely, that the spallation prod- 
ucts are not found initially in the accelerated CRs) are: 



i{t,t') = exp (-At/re )Ai 



(20) 



Appendix A: The "Age" of cosmic rays in the spiral 
arm diffusion model 



The tightest constraint on different types of CR dif- 
fusion models in the galaxy can be placed by compar- 
ing their predictions to the expected ratio between ^°Be, 
which is a radioactive isotope, to ^Be, which is stable. 
Both nucleotides are formed by CR spallation of mainly 
C, N and O. 

In a leaky box model, the ratio at steady state between 
the density of specie i of a stable isotope to that of an 
unstable isotope j, can be obtained from eq. ||. It is: 



Qi Te,lb + Tr 



(18) 



where Tr is the decay rate of the unstable isotope and 
Tejb is the leaky box decay rate. It is assumed that the 
spallation time is much longer than Te^w or r^. 

Different mea surem ents for the isotope ratio ^°Be/Be 
are 6.4 ± 1.5% |ll§, 3.9 ± 1.4% @, 4.3 ± 1.5% [||. 
These ratios respectively yield a Te^ib of 18.1^5' J Myr, 
30.3j:?J;^ Myr, and 27.1^^^ Myr ||. We shah adopt 



an age of 24^g^ Myr which is the measurements' average 
and a somewhat conservative error which is the maximum 
error divided by \/3. 
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(21 

with T^g = ^ +T^ ^, and At = t — t'. If we write eq. [T 
as: 



<^{x,z = 0,t)= / dt'G{x,t,t')d4{0,exp[~T/{At + St)]) 

J —00 

(22) 

with 5t = all {2D), then: 

N.^qif dt'G{x,t,t')exp[-{t-t' + 6t)/Te] (23) 
J —00 

and 

Nj « Qj [ dt'G{x,t,t') 
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At + St\ ( At + 6t 
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V Toff 



(24) 



If we define TZ = Nj/Ni, and compare the above result 
to those of the leaky box model, then the effective leaky 
box decay time Te^ib is related to Te and TZ through: 



Te.lb.eff ~ ^'^^ ^ 



(25) 
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This allows us to calculate the efTective age of the CRs 
that the Be isotope ratio would correspond to in a "leaky 
box" model. 



Appendix B: Systematic errors from diffusion and 
epicycles in the solar orbit 

In this work, we compare two completely different 
types of measurements of the spiral pattern speed. On 
one hand, there are the present day astronomical mea- 
surements of the spiral pattern speed. On the other hand, 
there are the temporal records of passages through spiral 
arms — either as recorded in the ice age reoccurrence, or 
as recorded in the variable CRF in Iron meteorites. The 
two types of measurements should yield the same result 
if the solar orbit is stable at its current radius. This is 
not the general case however, we are therefore required 
to estimate the error that can arise when comparing the 
two. 

Although collisions that drastically change the solar 
orbit are statistically unexpected, the Solar system un- 
dergoes many small angle deflection such that its orbit 
parameters slowly diffuse. Namely, the "Lagrangian" 
record of spiral arm crossings, as measured by the mov- 
ing solar system, need not be exactly the same as a test 
"Eulerian" particle would experience when fixed at the 
current Rq. 

The solar motion itself can be broken down to a cir- 
cular motion of the guiding center at radius Rm and an 
epicyclic motion at frequency k and radius an. 



with 



According to Wielen and Wielen et al. [|112[ , both 
amplitudes undergo diffusion, according to: 



Ai?n 



and 



= ^{ARD = ■ An = (0.72 kpc) VGyr, 

(26) 



{Aa\) = VDJ ; Da = (0.84 kpc)VGyr, 

(27) 

for t > 10^ yrs. Also relevant, is the "diffusion" along 
the azimuthal angle direction. This is given by: 



A5„ 



(12.7 kpc)VGyr3. 

(28) 

If we look at a given time t before present, the uncer- 
tainties in Rm and S will be given by the above formulae. 
The uncertainty in a_R is a trifle more complicated. This 
is because we have the prior knowledge that at — A.h 
Gyr BP, a^^Q = while today it is a^^i = 0.6 kpc. To 
calculate AaR.rmsj we calculate the probability of the so- 
lar system having a given a^^ at time t before present. It 
is given by: 



(29) 



G(aR, v/Da(to - t))G(aR. - aR.,i, - 
J^ao G(aR, VDa(to - t))G(aR - aR,i, 



G(x,(t) 



Using this, we find that: 



(31) 



J — OG 

l{to-t){Datot+{to-t)al^ 



This result is the expected amplitude of the epicyclic mo- 
tion as a function of time before present. 

The diffusion processes will translate into two types of 
errors. One is the measurement error in the location of 
each spiral arm. The second is an error arising when cal- 
culating the average pattern speed over several periods. 

The main contribution to the error in each spiral arm 
location is due to the epicyclic motion of the solar sys- 
tem around the 'guiding center', at a typical period of 
200 Myrs. This implies that each spiral arm crossing is 
inaccurately predicted. This translates into a phase shift 
within the spiral passages. The r.m.s. of the variation is: 



A(/)a 



AaR,rms/\/2 

ttRq sin i ' 



(32) 



assuming 4 arms inclined at an angle i, and that the 
r.m.s. of the "circular" epicyclic motions is l/\/2 of the 
radius. For a 140 Myr period (which will be a typical 
value found) , we obtain that the phase error ranges from 
A0nns = 15% today to A0„„s = 22% at 800 Myr BP. 
For a typical spiral arm crossing period of 140 Myr, this 
translates into an error of 21 to 31 Myrs respectively. 

If we measure the pattern speed by dividing the dura- 
tion it takes for TV spiral passages, several errors appear. 
First, the solar orbit can diffuse in or out. This intro- 
duces two errors by itself. The first is a phase error. 
Averaged over N spiral passages, it introduces a relative 
error of: 
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(33) 



where AfJin^rms is the value expected after N spiral pas- 
sages. This will give an absolute error of: 



^*/*m,rms, 1 







rip 



ttNRq tan I 
0.4 km sec^^ kpc^^ 



for t = 800 Myr 



in the relative pattern speed. The second error from the 
diffusion of the solar orbit arises from the fact that at 
a different orbit, the pattern speed it different. Again, 
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averaged over N spiral passages, it introduces an absolute 
error in the pattern velocity of: 



AR 



in,rms^0 

2R 

1.1 km scc"^ kpc~^ for t = 800 Mj^8.4) 



A phase error also arises from the diffusion along the 
orbit: 



Ar^s^rms — A0s,rms,l l^^i 



ttNRq 

« 0.3 km sec"^ kpc"^ for t = 800 Myr.(35) 

The first two errors are maximally correlated (they would 
have been anti correlated if either Qp > VIq or i < 0, but 
not both), and they are almost uncorrelated with the 
last. Therefore, the total error is: 



1.5 km sec"^ kpc'^ for t = 800 Myr.(36) 



Some indication to the radial solar migration can be 
obtained from its anomalously higher metallicity than its 
neighbors. This is statistically not very significant, but it 
does bias our expectation for solar diffusion. By compar- 
ing the solar metallicity to the surroundings and taking 
into account the apparent rad ial g radient in metallicity in 
the Milky Way, Wielen et al. ]112| found that the sun has 
migrated outwards by ARq = 1.9 ± 0.9 kpc. Namely, we 
expect that the probability for the formation radius Rq of 
the sun be given by G(R0 — Ro — 1.9 kpc, ctq — 0.9 kpc). 
Therefore, the probability at time t before present, to 
find the solar system at a radius R is given by: 



P{R,t) 



C^K{Ro,R,RQ,t)dRo 
J^^KiRo,R,RQ,t)dRodR 



(37) 



K{Ro,R,RQ,t) = G(R0-Ro-ARo,<To) 



G(R-Ro, VDm(to-t)) 
xG(R0 -R,v/D;;;t) (38) 



We then obtain that the expectation value for R at time 
t. It is: 



(R) 



Re 



RP{R,t)dR 

DmtARo 
Dmto + CTg 



(39) 



=0.8 Gyr Rq - 0.25 kpC 



and its r.m.s. error at t = 800 Myr is: 



Ai??^, ^ - {Rf 



(40) 



R^P{R,t)dR 



RP{R,t)dR 



D^t 1 



t=0.8 Gyr (0.59 kpc)^ 

Using the previous results, we obtain that for a mea- 
surement baseline spanning 800 Myrs BP until today, 
the difference between the "Lagrangian" and "Eulerian" 
measurements is 



(ri0-ilp)Lag-(ri0-rip)Euior ~ 0.54±1.5km s ^ kpc ^ 

(41) 



[1] Alvarez, L. W., Alvarez, W., Asaro, F., & Michel, H. V. 

1980. Science, 208, 1095. 
[2] Amaral, L. H., & Lepine, J. R. D. 1997. 

Mon. Not. Roy. Astr. Soc, 286, 885. 
[3] Avedisova, V. S. 1989. Astrophys., 30, 83. 
[4] Barry, D. C. 1988. Astrophys. J. , 334, 436. 
[5] Bazilevskaya, G. A. 2000. Sp. Sci. Rev., 94, 25. 
[6] Beer, J., Mende, W., & Stellmacher, R. 2000. Quat. Sci. 

Rev., 19, 403. 

[7] Begelman, M. C, & Rees, M. J. 1976. Nature (London) 
, 261, 298. 

[8] Bell, A. R. 1978. Mon. Not. Roy. Astr. Soc, 182, 443. 
[9] Berezinskii, V. S., Bulanov, S. V., Dogiel, V. A., & 
Ptuskin, V. S. 1990. Astrophysics of cosmic rays. Ams- 
terdam: North-Holland, 1990, edited by Ginzburg, V.L. 
[10] Binney, J., & Tremaine, S. 1988. Galactic Dynamics. 

Princeton University Press, Princeton. 
[11] Blitz, L., Fich, M., & Kulkarni, S. 1983. Science, 220, 
1233. 



[12] Christensen-Dalsgaard, J., Dilke, F. W. W., & Gough, 
D. O. 1974. Mon. Not. Roy Astr. Soc, 169, 429. 

[13] Comeron, F., & Torra, J. 1991. Astron. Astrophy, 241, 
57. 

[14] Compton, A. H. 1933. Phys. Rev., 43, 387. 
[15] Creze, M., & Mennessier, M. O. 1973. Astron. Astro- 
phy, 27, 281. 

[16] Crowell, J. C. 1999. Pre-Mesozoic Ice Ages: Their Bear- 
ing on Understanding the Climate System. Vol. 192. 
Memoir Geological Society of America. 

[17] Dame, T. M., Hartmann, D., & Thaddeus, P. 2001. As- 
trophys. J. , 547, 792. 

[18] Dickinson, R. E. 1975. Bull. Amer. Met. Soc, 56, 1240. 

[19] Dilke, F. W. W., & Gough, D. O. 1972. Nature (Lon- 
don) , 240, 262. 

[20] Dopita, M. A. et al. 1997. Astrophys. J. , 474, 188. 

[21] Dragicevich, P. M., Blair, D. G., & Burman, R. R. 1999. 
Mon. Not. Roy Astr. Soc, 302, 693. 

[22] Duric, N. 2000. Supernova remnants and cosmic rays 



32 



in M31 and M33. Pages 179-186 of: Proceedings 232. 
WE-Heraeus Seminar, 22-25 May 2000, Bad Honnef, [55; 
Germany. Edited by EUy M. Berkhuijsen, Rainer Beck, [56 
and Rene A. M. Walterbos. Shaker, Aachen, 2000, p. 
179. [57 
Efremov, Y. N. 1983. Sov. Astron. Lett., 9, 51. 
Egorova, L. Y., Vovk, V. Ya., & Troshichev, O. A. 2000. 
J. Atmos. Solar-Terr. Phys., 62, 955. 
Elmegreen, D. M. 1998. Galaxies and Galactic Struc- 
ture. Prentice Hall. [58 
Evans, N. J. 1991. in ASP Conf. Ser. 20: Frontiers of [59; 
Stellar Evolution, p. 45. [60^ 
Fastrup, B. et al. 2001. A study of the link between 
cosmic rays and clouds with a cloud chamber at the [61 
CERN PS. CERN-SPSC-2000-021. Also as: LANL 
physics/0104048. [62; 
Feitzinger, J. V., & Schwerdtfeger, H. 1982. Astron. As- 
trophy., 116, 117. 
Fernandez, D., Figueras, F., & Torra, J. 2001. As- [63 
tron. Astrophy., 372, 833. 

Frakes, L. A., Francis, E., & Syktus, J. I. 1992. Climate [64 
modes of the Phanerozoic; the history of the Earth's 
climate over the past 600 million years. Cambridge: 
Cambridge University Press, 1992. [65 
Friis-Christensen, E., & Lassen, K. 1991. Science, 254, 
698. [66 
Gallagher, J. S. et al. 1996. Astrophys. J. , 466, 732. [67 
Gardiner, L. T., Sawa, T., & Fujimoto, M. 1994. 
Mon. Not. Roy. Astr. Soc, 266, 567. [68; 
Georgelin, Y. M., & Georgehn, Y. P. 1976. Astron. As- 
trophy., 49, 57. [69 
Gordon, M. A. 1978. Astrophys. J. , 222, 100. 
Grivucv, E. M. 1981. Sov. Astron. Lett., 7, 303. [70' 
Grivucv, E. M. 1983. Sov. Astron. Lett., 9, 287. 
Guyodo, Y., & Valet, J. 1996. Earth Planet. Sci. Lett., [71 
143, 23. 

Haigh, J. D. 1996. Science, 272, 981. [72; 

Hambrey, M. J., & Harland, W. B. 1985. Palaeogcog- 

raphy, Palaeoclimatorlogy, Palacoocology, 51, 255. [73 

Hampel, W., & Schacffer, O. A. 1979. Earth. Planet. 

Sci. Lett., 42, 348. [74 

Harrison, R. G. 2000. Sp. Sci. Rev., 94, 381. [75^ 

Hernandez, X., Valls-Gabaud, D., & Gilmore, G. 2000. 

Mon. Not. Roy. Astr. Soc, 316, 605. [76' 

Herschcl, W. 1796. Phil. Trans. Roy. Soc, London, 166. 

Herschel, W. 1801. Phil. Trans. Roy. Soc, London, Part [77 

1, 265. 

Heyer, M. H., Brunt, C, Snell, R. L., Howe, J. E., Schlo- [78 
erb, F. P., & Carpenter, J. M. 1998. Astrophys. J. Supp., [79^ 
115, 241. 

Hodell, D. A., Brenner, M., Curtis, J. H., & Guilderson, [80' 
T. 2001. Science, 292, 1367. [81 
Hoffman, P. F., Kaufman, A. J., Halverson, G. P., & 
Schrag, D. P. 1995. Science, 281, 1342. [82; 
Hoyle, F., & Lyttleton, R. A. 1939. Proc. Cambridge 
Phil. Soc, 35, 405. [83; 
Hoyle, F., & Wickramasinghe, C. 1978. Astro- 
phys. Space Sci., 53, 523. [84 
Huber, B. T. 1998. Science, 282, 2199. 
IPCC. 1995. Climate Change 1995, Intergovernmental 
Panel on Climate Change. Cambridge University Press. [85 
Ivanov, G. R. 1983. Pis ma Astronomicheskii Zhurnal, 
9(Apr.), 200. [ 
Kirkby, J., & Laaksonen, A. 2000. Sp. Sci. Rev., 94, 



397. 

Labitzke, K., & van Loon, H. 1992. J. Glim., 5, 240. 
Lavielle, B., Marti, K., Joannot, J., Nishiizumi, K., & 
Caffee, M. 1999. Earth Planet. Sci. Lett., 170, 93. 
Lees, J. F., & Lo, K. Y. 1990 (July). Dissociation and 
ionization of molecular gas in the spiral arms of M51. in 
NASA, Ames Research Center, The Interstellar Medium 
in External Galaxies: Summaries of Contributed Pa- 
pers, p 296 (SEE N91-14100 05-90). 
Leitch, E. M., & Vasisht, G. 1998. New Astron., 3, 51. 
Leitherer, et al. 1999. Astrophys. J. Supp., 123, 3. 
Lin, D. N. C, Jones, B. F., & Klemola, A. R. 1995. 
Astrophys. J. , 439, 652. 

Lisonfcld, U., Alexander, P., Pooley, G. G., & Wilding, 

T. 1996. Mon. Not. Roy. Astr. Soc, 281, 301. 

Lo, K. Y., Ball, R., Masson, C. R., Phillips, T. G., Scott, 

S., & Woody, D. P. 1987. Astrophys. J. Letters, 317, 

L63. 

Loktin, A. V., Matkin, N. V., & Gerasimenko, T. P. 
1994. Astron. Astrophys. Trans., 4, 153. 
Longair, M. S. 1994. High energy astrophysics. Vol.2: 
Stars, the galaxy and the interstellar medium. Cam- 
bridge: Cambridge University Press, 1994, 2nd ed. 
Lukasiak, A., Ferrando, P., McDonald, F. B., & Webber, 
W. R. 1994. Astrophys. J. , 423, 426. 
Lynds, B. T. 1980. Astron. J., 85, 1046. 
Marsh, N., & Svensmark, H. 2000. Sp. Sci. Rev., 94, 
215. 

Michel, R. et al. 1995. Nuc. Inst. Meth. Phys. Res. B, 
103, 183. 

Mishurov, I. N., Pavlovskaia, E. D., & Suchkov, A. A. 
1979. Astronomicheskii Zhurnal, 56, 268. 
Mishurov, Y. N., & Zenina, I. A. 1999. Astron. Astro- 
phy., 341, 81. 

Napier, W. M., & Clube, S. V. M. 1979. Nature (Lon- 
don) , 282, 455. 

NefF, U., Burns, S. J., Mangnini, A., Mudelsee, M., 

Fleitmann, D., & Matter, A. 2001. Nature, 411, 290. 

Nelson, A. H., & Matsuda, T. 1977. 

Mon. Not. Roy. Astr. Soc, 179, 663. 

Ney, E. P. 1959. Nature, 183, 451. 

Olhng, R. P., & Merrifield, M. R. 1998. 

Mon. Not. Roy. Astr. Soc, 297, 943. 

Palle Bago, E., & Butler, J. 2000. Astron. & Geophys., 

41, 18. 

Palous, J., Ruprecht, J., Dluzhnevskaia, O. B., & 
Piskunov, T. 1977. Astron. Astrophy., 61, 27. 
Perko, J. S. 1987. Astron. Astrophy., 184, 119. 
Pudovkin, M. I., & Veretenenko, S. V. 1995. J Atmos. 
Terr. Phys., 57, 1349. 

Rind, D., & Overpeck, J. 1993. Quat. Sci. Rev., 12, 357. 
Robin, A. C, Creze, M., & Mohan, V. 1992. Astro- 
phys. J. Letters, 400, L25. 

Rocha-Pinto, H. J., Scalo, J., Maciel, W. J., & Flynn, 
C. 2000a. Astron. Astrophy., 358, 869. 
Rocha-Pinto, H. J., Scalo, J., Maciel, W. J., & Flynn, 
C. 2000b. Astrophys. J. Letters, 531, L115. 
Ruphy, S., Robin, A. C, Epchtcin, N., Copet, E., 
Bertin, E., Fouque, P., & Gugliclmo, F. 1996. As- 
tron. Astrophy., 313, L21. 

Scalo, J. M. 1987. in Starbursts and Galaxy Evolution, 
p. 445. 

Schaeffer, O. A., Nagel, K., Fechtig, H., & Neukum, G. 
1981. Planet. Space Sci., 29, 1109. 



33 



[87] Shaviv, N. J. 2002. Phys. Rev. Lett. , 89, 051102. 

[88] Simpson, J. A., & Garcia-Munoz, M. 1988. Sp. Sci. 

Rev., 46, 205. 
[89] Singer, S. F. 1954. Astrophys. J. , 119, 291. 
[90] Soon, W. H., Posmentier, E. S., & Baliunas, S. L. 1996. 

Astrophys. ,J. , 472, 891. 
[91] Soon, W. H., Posmentier, E. S., & Baliunas, S. L. 2000. 

Annalcs Geophysicao, 18, 583. 
[92] Steiner, J., & Grillmair, E. 1973. Gcol. See. Am. Bull., 

84, 1003. 

[93] Stozhkov, Yu. I., & et al. 1995. II Nuovo Cimento C, 
18, 335. 

[94] Svensmark, H. 1998. Phys. Rev. Lett. , 81, 5027. 

[95] Svensmark, H., & Friis-Christensen, E. 1997. J. Atmos. 
Terr. Phys., 59, 1225. 

[96] Tammann, G. A., Loffler, W., & Schroder, A. 1994. As- 
trophys. J. Supp., 92, 487. 

[97] Taylor, J. H., & Cordes, J. M. 1993. Astrophys. J. , 411, 
674. 

[98] Tinsley, B. A., & Been, G. W. 1991. J. Geophys. Res., 

12, 22,283. 

[99] Vallec, J. P. 1995. Astrophys. J. , 454, 119. 
[100] Vallee, J. P. 2002. Astrophys. J. , 566, 261. 
[101] Vallenari, A., Chiosi, C., Bertelli, G., & Ortolani, S. 

1996. Astron. Astrophy., 309, 358. 
[102] van den Bergh, S., & McClure, R. D. 1994. Astrophys. 

J. , 425, 205. 

[103] Veizer, J., Godderis, Y., & Francois, L. M. 2000. Nature, 
408, 698. 

[104] Voshage, H. 1967. Zeitschrift Naturforschung Teil A, 
22, 477. 

[105] Voshage, H., & Feldmann, H. 1979. Earth Planet. Sci. 



Lett., 45, 293. 

[106] Voshage, H., Fcldmanu, H., & Braun, O. 1983. 

Zeitschrift Naturforschung Teil A, 38, 273. 
[107] Wagner, G., Beer, J., Kubik, P. W., & Synal, H.-A. 

1999. in Annul Report of the Paul Sclicrrcr Institute, 

Swiss Federal Institute of Technology (ETH) , Zurich. 
[108] Webber, W. R., & Soutoul, A. 1998. Astrophys. J. , 

506, 335. 

[109] Wosterlund, B. E. 1990. Astron. Astrophys. Rev., 2, 
29. 

[110] Wiedenbeck, M. E., & Greiner, D. E. 1980. Astro- 
phys. J. Letters, 239, L139. 

[Ill] Wielen, R. 1977. Astron. Astrophy., 60, 263. 

[112] Wielen, R., Fuchs, B., & Dettbarn, C. 1996. Astron. As- 
trophy., 314, 438. 

[113] Williams, G. E. 1975. Earth Planet. Sci. Lett., 26, 361. 

[114] Yabushita, S., & Allen, A. J. 1985. The Observatory, 
105, 198. 

[115] Yu, F. 2002. J. Geophys. Res., in press, . 

[116] Yuan, C. 1969a. Astrophys. J. , 158, 871. 

[117] Yuan, C. 1969b. Astrophys. J. , 158, 889. 

[118] Zachos, J., Pagani, M., Sloan, L., Thomas, E., & 
Billups, K. 2001. Science, 292, 686. 

[119] A common misconception is that the spirals are "frozen" 
in, such that material in the arms remains in them, and 
vice versa. If this would have been true, the spirals of 
galaxies should have boon much more tightly wound be- 
cause of the differential rotation. 

[120] Because of limited coverage in Galactic longitude, 3 spi- 
ral arms are seen. 2 adjacent arms end within the cov- 
ered longitude, at the same inferred radius of i?Hi,out 'a 
2Rq . The separation angles imply a ~ 90° separation. 



